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Abstract 

Hadron spin physics is now one of the most active fields of physics. Especially in 
the last ten years, great progress has been made both theoretically and experimentally 
that has considerably improved our knowledge of the spin structure of nucleons. We 
review the nucleon's polarized structure functions from the viewpoint of factorization 
theorems and the gauge invariant, nonlocal light-cone operators in QCD. We discuss a 
systematic treatment of the polarized structure functions and the corresponding parton 
distribution functions, which are relevant to inclusive lepto- and hadro-production. We 
give a detailed analysis of these spin-dependent distribution functions at the twist-2 
and twist-3 level, and present various properties and relations satisfied by the parton 
distributions, which can be derived directly from QCD. We emphasize unique features 
of higher twist distributions, and the role of the QCD equations of motion to derive 
their sensitivity to the quark-gluon correlation and their anomalous dimensions for 
Q 2 -evolution. 
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1. Introduction 



Hadron spin physics is now one of the major fields of particle and nuclear physics based 
on quantum chromodynamics (QCD). The starting point for this field was the measurement 
of the polarized nucleon structure function g\(x,Q 2 ) in deep inelastic scattering by the 
EMC collaboration. & The data of this EMC experiment indicate that an unexpectedly 
small fraction of the nucleon's spin was carried by its constituents, quarks, and have caused 
many physicists to challenge the so-called "spin crisis" problem in QCD. After a flood of 
theoretical papers as well as experimental data,B ) 'i''B ) our understanding on this problem 
is much more updated. 0* Through these developments, "hadron spin physics" has grown up 
as a field attracting considerable attention. Now our interest has spread to various other 
processes to explore the spin structure of the nucleon. In particular, in conjunction with new 
projects like the "RHIC spin project", "polarized HERA", etc., we are now in a position to 
obtain more information on the deep structure of the nucleon and the dynamics of QCD. 
Furthermore, the spin degrees of freedom in high-energy processes are expected to play an 
important role to verify the Standard Model itself. 

The aim of this article is to provide an overview of the perturbative methods in QCD 
used to study the nucleon's spin structure through deep inelastic processes, and to summa- 
rize the recent theoretical progress on this subject including the QCD evolutions of polarized 
structure functions. Although the high-energy (short- distance) behavior of QCD is de- 
scribed by perturbation theory thanks to asymptotic freedom, a general cross section is a 
complicated combination of the short- and long-distance interactions, and is not directly cal- 
culable. For some specific processes, however, the factorization theorems allow us to separate 
(factorize) the short-distance dynamics from the long-distance dynamics in the cross section 
in a systematic fashion, and to derive the predictions for the short-distance dynamics using 
perturbation theory. Many observables are known in spin related phenomena of high-energy 
processes, to which factorization theorems and thus perturbative QCD can be successfully 
applied. In this paper, we restrict ourselves to the nucleon's polarized structure functions, 
which can be measured in high-energy inclusive reactions, like lepton-nucleon deep inelastic 
scattering and the Drell-Yan process in nucleon-nucleon collisions. 

In the present context, the proper high-energy limit is provided by the "Bjorken limit" . In 
this limit, using the notion of "twist" , we can systematically expand a cross section in inverse 
powers of a characteristic large momentum scale Q, and extract the leading contributions. 
The factorization theorems can be proved (in principle) order by order in the twist expansion. 

*> Those who are interested in the present status of phenomenological studies are recommended to refer 
to the recent comprehensive review by Lampe and Reya. lk 
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At the leading twist (twist-2) level, the result implied by the factorization theorems is realized 
as the QCD improved parton model. On the other hand, the higher twist effects describe 
the coherent quark-gluon behavior inside the nucleon, and thus contain information on the 
correlation of quarks and gluons beyond the parton model. The higher twist contributions 
are suppressed by some powers oil/Q compared to the leading twist contribution and are 
usually hidden by the latter in the cross section. This makes it very difficult to extract 
higher twist effects from experimental studies. In the case of spin-dependent processes, 
however, we have a good chance to study them; it is known that the twist-3 contributions 
can be measured as leading effects in certain asymmetries. Therefore we will discuss in detail 
the twist-3 contributions in addition to the leading twist effects for the polarized structure 
functions. 

This paper is organized as follows. In §2 we recall the factorization theorems in the 
above two typical examples of inclusive processes. To describe the long-distance part of 
the factorized cross section, we introduce the parton distribution functions as the nucleon's 
matrix elements of nonlocal light-cone operators in QCD. We discuss the parton distributions 
given by the two- and three-particle operators and classify them with respect to twist, spin- 
dependence and chiral properties. In §3 we explain the equivalence of the approach in §2 
to the conventional approaches based on the operator product expansion. A general feature 
of the renormalization of gauge invariant operators will be discussed in §4 to deal with the 
higher twist operators in a covariant approach. The QCD predictions for the chiral-even 
structure functions are summarized in §]S|. The chiral-odd structure functions are discussed 
in §EL The final section is devoted to a summary. 



§2. Factorization theorems and parton distribution functions 

To describe a variety of high-energy processes in a universal language, it is desirable 
to have a definition of parton distribution functions based on the operators in QCD. The 
traditional approach for this purpose relies on the operator product expansion (OPE), but it 
can be applied only to a limited class of processes like deep inelastic lepton-hadron scattering. 
This calls for an approach based on the factorization as a generalization of the OPE. In this 
section, we explain the factorization theorems, and introduce a definition of the parton 
distribution functions in terms of the nonlocal light-cone operators in QCD. 

2.1. Factorization theorems for hard processes in QCD 

There are various hard processes that are characterized by the large momentum squared 
Q 2 : deep inelastic scattering (DIS), l+A —> I'+X (A, B, . . . represent hadrons, I a lepton, and 
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X a system of hadrons produced through inelastic processes); the Drell-Yan (DY) processes, 
A + B — > l + + P + X ; jet production, A + B — > jet + X ; heavy quark production, 
A + B — > heavy quark + X, etc. The basis for the QCD analysis of these hard processes is 
provided by the factorization theorems in QCD,& which give an extension of the OPE. It 
provides a foundation of the "parton model" in the Bjorken limit, where Q 2 — > oo with the 
Bjorken variable x fixed, and also a systematic framework to calculate the QCD corrections 
beyond the leading order. The theorems tell us that we can view a high-energy beam 
of hadrons as if it were a beam of partons (quarks and gluons), and the hadron reaction 
is induced by hard scattering among individual partons. Corresponding to this intuitive 
picture, the cross section for the hard processes is given as a product, or more precisely, 
convolution of the short- and long-distance parts. The former contains all the dependence 
on the large momentum Q, while the latter depends essentially on the QCD scale parameter 
Aqcb- These two parts are divided at a factorization (renormalization) scale and the 
short-distance (long-distance) part involves the momenta larger (smaller) than \i. The short- 
distance part corresponds to the hard scattering cross section for the partons with the large 
momentum Q exchanged, while the long-distance part corresponds to the parton distribution 
functions in a hadron. The former is systematically calculable in perturbation theory for 
each process, due to the asymptotic freedom of QCD, while the latter is controlled by the 
nonperturbative dynamics of QCD. 

Let us recall the factorization formulae for some familiar processes. First, consider the 
DIS, which proceeds via the exchange of a virtual photon with momentum between a 
lepton and a hadron (Q 2 = —q 2 > 0). S As is well known, the DIS cross section is given in 
terms of the leptonic (L^) and the hadronic (W^) tensors, H* 



da _ 1 / e 2 \\^ w 

where k(k') and P denote the momenta of the incident (scattered) lepton and hadron (k—k' 
q) . All the information of the strong interaction is contained in the hadronic tensor: 

WW = 7T T,( PS W)\X)(X\JM\PS)(2nr5^\p x -P-q) 
2tt x 



1 

2^ 



d A ze-^ z {PS\[J^) , Jv{z))\PS) . (2-1) 



*> The factorization scale can be different from the renormalization scale. In this paper, we set them 
equal for simplicity. 



**) More generally, the vector boson or Z° can be exchanged. In this case, J M (z) of (2-2) should be 
replaced by the corresponding charged or neutral current. 
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Here \PS) is the hadron state with momentum P and spin S, px is the momentum of the 
hadronic final state \X), and 

J,(z) = $(z){\ 3 /2 + A 8 /(2v / 3)} 7 ^(^) = Hz)& el hMz) (2-2) 

is the hadron's electromagnetic current composed of the quark field ip(z) an d the flavor 
matrices A* for the u-, d- and s-quarks when we consider three flavors. For the case of a 
nucleon target with mass M (P 2 = M 2 , P ■ S = 0, S 2 = —M 2 ), the hadronic tensor can be 
expressed in general by the four structure functions Fi,F 2 , g\ and g 2 as 

W vv = 

- ( ^ - q -f) Q 2 ) + ( p> - ( p» - ^) j?r q Wx, g 2 ) 

+ 2t e, uXa q x j^p^ 9i(x, Q 2 ) + (P " QS a - q ■ S P°) -^-^ g 2 (x , Q 2 ) | . (2-3) 

The symmetric (antisymmetric) part in \iv is relevant to the unpolarized (polarized) scat- 
tering; in (2-3), the decomposition of the spin-independent part is universal for an arbitrary 
hadronic target, while that of the spin-dependent part is specific for the spin 1/2 target. 
The structure functions are dimensionless functions of two invariants, Q 2 and the Bjorken 
variable x = Q 2 /(2P ■ q). Equation (2-1) can be expressed by the "cut diagram" corre- 
sponding to the discontinuity of the forward virtual Compton amplitude between the virtual 
photon and a hadron. In general kinematics, this virtual Compton amplitude contains all 
the complicated interactions between the virtual photon and a hadron, including various 
"soft interactions" where soft momenta are exchanged. However, a drastic simplification 
occurs if one goes to the Bjorken limit Q 2 — > oo with x fixed: the amplitude is dominated 
by the contribution which is factorized into the short- and long-distance parts (see Fig. 1), 
and other complicated contributions are suppressed by the powers of 1/Q. 




Fig. 1. The DIS process and the Bjorken limit. 

The factorized amplitude corresponds to the process in which a parton carrying the 
momentum £P (0 < £ < 1) comes out from the long-distance part, followed by the hard 
impact of the virtual photon, and then goes back to the long-distance part. As a result, the 
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structure functions F±(x, Q 2 ) and F 2 (x, Q 2 ), which give the unpolarized DIS cross section up 
to kinematical factors, assumes the factorized form 



x q_ 



(24) 



and similarly for Fi(x,Q 2 ). The quantity fi/Ai^fJ 2 ) is the parton distribution function 
corresponding to the long-distance part, and it is interpreted as the probability density of 
finding a parton of type i (= gluon, u, u, d,d,- ■ •) in a hadron A, carrying a fraction £ of the 
hadron's momentum. The summation of ( |24| ) is over all possible types of partons, i. Hi 
denotes the short- distance part; it corresponds to the hard scattering cross section between 
the virtual photon and a parton i, and is given by a power series in a s = g 2 j^n (g is the 
QCD coupling constant) with finite coefficients. 

Next, we consider the DY process. & The cross section for this process is again written 
in terms of the leptonic and hadronic tensor, 



da 



s(s-AM 2 ) \Q 



JiY r^w d3kl d3ft2 

-2 j *- Kl V (2^)3 2^ (2tt) 3 2A; 2 ' 



where ki is the momentum of the produced lepton, Q 2 = q 2 = [ki + k 2 ) 2 and s = (Pa + Pb) 2 - 
The hadronic tensor is given byi'i) 



VIV = J d A ze~^ z (P A SA,PBSB\J^z)JM\PASAiPBS B ) 



(2-5) 



where \PaSa, PbSb) is the in-state of two hadrons A and B. g^ u W^ v determines the DY 
cross section with unpolarized hadron-hadron collision, while the other components of W^ u 
are relevant to polarized collisions. Equation (|2-5|) is expressed by the cut of the forward 
scattering amplitude between A and B (see Fig. 2). In the Bjorken limit Q 2 — ► oo with 
r = xaXb = Q 2 /s fixed, the amplitude is dominated by the contribution factorized into 
short- and long-distance parts; other contributions are suppressed by the powers of l/Q. 




Q 2 ^oo 

x : fixed 




Fig. 2. The DY process and the Bjorken limit. 
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The factorized amplitude corresponds to the process in which partons carrying the momenta 
£aPa and £bPb come out from the lower and the upper long-distance parts, followed by 
the hard scattering between them, and then go back to the long-distance parts (see Fig. 2). 
Then, the cross section for the unpolarized DY, accurate up to corrections suppressed by 
the powers of l/Q, is given by 



d 2 



a 



dQHy 



I [ d£ B f i/A (£ A , fi)H i:i I ^-,Q, — ,a s (/i 2 ) J fj/ B (^B,fj) 



(2-6) 



where xa/xb = e 2y . Here is the hard-scattering coefficient corresponding to the short- 
distance interaction between the partons i and j. The long-distance parts appearing in Fig. 2 
have exactly the same structure as that appearing in Fig. 1, and therefore fa/ a and fj/B are 
the same as in (2_4). This demonstrates the universal nature of the long-distance part: it 
is determined completely if one specifies the target. Though the long-distance part is not 
calculable using the perturbation theory, it can be determined by experiments of some hard 
processes. For example, if one extracts the parton distribution functions for the nucleon 
target by DIS experiments, they can be used to describe the DY processes involving the 
nucleon. 

The validity of ( |2-4j ) and ( |2-6| ) can be demonstrated by analyzing the Feynman diagrams 
for the corresponding amplitude in the Bjorken limit. & (For the DIS, the equivalent result 
can be obtained by the OPE. I) The central point in the proof is to show that the long- 
distance contributions, which appear as infrared- divergent contributions in the perturbative 
diagrams, can be completely factorized into universal hadron matrix elements (parton dis- 
tribution functions) or cancel out, and thus the hard scattering coefficient has purely short- 
distance contributions. Extensions of the factorization theorems to hard processes involving 
polarized beams and/or targets have been discussed by many authors. The correspond- 
ing factorization formulae are similar to ( p-4[) and ( |2-6|) but involve various polarized parton 
distribution functions. 



2.2. Factorization of the DIS in the free field theory 

Once the factorization theorems are proved, the parton distribution functions can be 
studied separately from the short-distance part and deserve detailed discussions because of 
their universal nature. Before going into a systematic study of parton distribution functions, 
it is instructive to demonstrate the DIS formulae like ( |2-4| ) in free field theories, that is, 
neglecting the QCD interaction. The result provides a guide to operator definitions of the 
parton distributions. 

Let us first recall some basic points, which generally hold even in the presence of QCD 
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interactions. Here and below, we consider the nucleon's structure functions, and \PS) de- 
notes the nucleon state. As mentioned above, the relevant quantity for the DIS through 
one photon exchange is the hadronic tensor (2-1). In the Bjorken limit, the hadronic tensor 
(2-1) is governed by the behavior of the current products near the light-cone. This light-cone 
dominance is easily understood from the following observation. We take the Lorentz frame 
in which the momenta of the nucleon and the virtual photon are both in the e z direction, 
without loss of generality. It is convenient to introduce two auxiliary light-like vectors, 

^ = -^(1,0,0,1), «^=* (1,0,0,-1), 

with p 2 = w 2 = , p ■ w = 1. We have P^ = pf + ^-u> M , and V = P + is a parameter which 
.elects a specihc frame in the direction of e, In the Bjorken limit, C? = oo with 

x = Q 2 /(2P ■ q) fixed, it is easy to derive 

lim g" = (p ■ q + ^M 2 x^j w» - xp? + O (-^j . 



Therefore, writing the coordinate as z^ = r/p 11 + Xw^ + Zj_, we obtain limBj q ■ z 
P ■ qt] — Ax. By the Riemann-Lebesgue theorem, only the integration regions \q ■ z\ £ 1, 
namely \r]\^l/(P ■ q) and |A| ^ contribute to (2-1). This is the light-cone dominance 
which implies that only the region 

0<z 2 = 2r / A-4<2|r / A|<^ (2-7) 

Q 

is important in the Bjorken limit. The lower limit of z 2 results from causality. 

Now we consider the free field theory. For simplicity, we neglect the flavor structure of 
quarks and also quark mass (we write the current as = ^j^i/j). The current commu- 
tator in (2-1) can be easily worked out by using the free- field anticommutation relation:!^ 

{V>(O),-0(»} = - /dA{z), where A{z) = -(l/2n)e(z°)5(z 2 ) with e(z°) = z°/\z°\. We 
obtain 

[J M (0) , J u (z)} = - (d a A(z)) [S^oUvfr z ) ~ ^ m uMl(0, z)\ , (2-8) 
with S^aua = g^aQva — Q^vQaa + fiVrfiW, and the "nonlocal (bilocal) operators" defined by 

Z4(0, z) = V(0) 7 ^W - ^(^7^(0) , (2-9) 
Ul(0, z) = ^(OhW 0) + Hz)^l 5 m . (2-10) 



*) The light-cone coordinates a M = (a + , a~, a±) are related to the usual coordinates as a* = (a°±a 3 )/v / 2, 
a± = (a 1 , a 2 ), so that a ■ b = a + b~ + a~b + — a± ■ b±. Note that {w^) has only p + (w~) as nonzero 
components in these coordinates. 

**) We keep only terms which contribute to the DIS, omitting the c-number contributions from the 
nonlocal operators U VA (0, z). Thus Uy A (0, z) should be understood to be normal ordered. 
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We take the nucleoli matrix element of ( |2-8|) and substitute it into (2-1). In accord with the 
general argument fl2-7|) , A(z) of ( |2-8| ) selects the integration region with rj ~ z± ~ in the 
Bjorken limit. Therefore, we can expand the matrix element in powers of the deviation from 
the light-cone, - Xw», and approximate it as (PS\U VA (0, z)\PS) ~ (PS\U VA (0, Xw)\PS) 
(the neglected higher-order terms give contributions suppressed by powers of 1/Q after the 
integration). 

Let us consider the contribution of Vfy(0, Xw) in detail: the matrix element corresponding 
to the first term of (2-9) can be parameterized as 

(PS\$(0)7 a ip(\w)\PS) = 2p a q(X) + 2w a M 2 f 4 {X) . (2-11) 

The second term has a factor M 2 to cancel the dimensionality of «A Here, let us use the 
term "twist" in the following sense: when a piece of the operator matrix element contributes 
to the hard processes at order (1/Q)' -2 , such a piece is said to have twist-t. Actually, this 
definition of the twist based on the power counting of 1/Q has a slight mismatch with the 
conventional definition as "dimension minus spin" of the relevant operators. However, 
it is useful to classify the operator contribution to the cross sections with respect to the 
power dependence in 1/Q. In the above example, the first term is the twist- 2 contribution, 
and the second term is the twist-4 contribution, since the factor M 2 should be compensated 
dimensionally as M 2 /Q 2 in the physical structure functions. 

We insert ( |2-8| ) and (|2T1| ) into (2-1), and compute the corresponding contribution to W^ v 
to leading order in 1/Q. By introducing the Fourier transformations 

A{z) = I i0y elk ' z [ 2 ™( k °) 5 ( k2 \ > 

q(X) = J d£e-^q(0 = J d^<q(0 , 
it is straightforward to obtain 

W, v = 2S mva f J d£(q + &) a 6(q Q + ip°)6({q + ^) 2 ) [q(£) - <?(-£)] + ■ • • , (2-12) 

where the dots stand for the contributions from /4(A) and ie IMava (PS\U A (0, Xw)\PS); those 
from /4(A) produce the contribution suppressed by M 2 /Q 2 , as noted above, and thus are 
irrelevant here. On the other hand, those from ie fJiaua (PS\U A (0, Xw)\PS) are relevant to the 
antisymmetric part of W^ u . Noting that (q + £p) 2 = q 2 + 2£p • q = 2P ■ q(£ - x) + 0(M 2 /Q 2 ), 
we obtain 



) The mismatch gives rise to the "Wandzura-Wilczek parts" for twist-3 contributions (see §§|| and 
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Comparing this result with the general expression (2-3), one obtains for the unpolarized 
structure functions, 

F^x) = -F 2 {x) = f ^5{x/i - 1) [q(0 - q(-0] = q(x) - q(-x) , (2-14) 
x J c, 

with 

q( x ) = J _e a ^(P»S|^(0) jh^{Xw)\PS) , (2-15) 

up to corrections of 0{M 2 /Q 2 ). We see that q(£) corresponds to a parton distribution 
/i/A(C)A* 2 ) °f (El)' while Hi(x/£,,Q 2 /[i 2 ,g(fi 2 )) — > 6(x/£ - 1) in the present free case. As 
discussed in §2.3 below, q(x) and — q{— x) give the quark and the anti-quark distribution 
functions, respectively. Similarly, it is straightforward to see that the terms corresponding 
to the dots in (|2T3|) give results for the polarized structure functions of (2-3) as 

gi (x) = -(Aq(x)+Aq(-x)), gi(x)+g 2 (x) = ^(g T (x)+g T (-x)) , (2-16) 

to 0(M 2 /Q 2 ) accuracy. The distribution functions Aq(x) and gr{x) are defined similarly to 
(|2 • 1 5|) , and their explicit forms are given in (2-19) and (2-23) below. 

Equations (|2-8| ) and ( |2-12j ) demonstrate that the singular function d a A(z) corresponds to 
the short-distance part containing all the dependence on g M , while the matrix element of the 
nonlocal operators corresponds to the long-distance part yielding parton distribution func- 
tions. Roughly speaking, the above results are modified in two respects in the presence of the 
QCD interaction. First, the higher order interactions produce the logarithmic (ln(Q 2 /// 2 )) 
corrections to A(z), and correspondingly the parton distribution functions acquire a depen- 
dence on the renormalization scale /i. Second, the coupling of the "longitudinal" gluons 
replaces A 

A(z)^A(z)[0,z] , 



where 



[y, z]=P exp {igj Q dt {y - z)^{ty + (1 - t)z)^j 



is the path-ordered gauge phase factor along the straight line connecting the points z^ and 
y^. By absorbing this factor, the nonlocal light-cone operators which define the parton 
distribution functions now preserve gauge invariance. 

Extension to 3 flavors with the current ( p^) is straightforward. The quark charge matrix 
squared, Q^ 2 = |i + |A 3 + A 8 , should be inserted between the quark fields of (2-9) 
- fl2Tl| ) and ( |2-15| ). The long-distance part is, therefore, decomposed into singlet (l) and 
non-singlet (A 3 ' 8 ) parts. 



The factor 1/2 on the r.h.s. comes from our definition of gi in (2-3). 
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2.3. Classification of the quark distribution functions 

Now, it is not difficult to extend the above logic in the free case to the quark distributions 
in the presence of QCD interaction.© We consider the following quantity involving the 
nucleon matrix element: 

1-°° dX 



2tt 



(p s \^(o)[o, Xw]r^(Xw)\PS) . 



(2-17) 



Here J 1 is a generic Dirac matrix, and the gauge phase factor [0, Xw] makes the operators 
gauge invariant. Equation fl2-17|) defines the distribution function for a quark with momen- 
tum k + = xP + (recall Ax = (xP) ■ {Xw)). It is understood that the nonlocal light-cone 
operator in (|2-17|) is renormalized at the scale \i. 

P can be any Dirac matrix, depending on which hard process is considered. An important 
observation made by Jaffe and Ji0) is that one can generate all quark distribution functions 
up to twist-4 by substituting all the possible Dirac matrices for r, although the operator 
definition (|2-17| ) is motivated by the factorization at the leading twist. By decomposing 
( |2-17p into independent tensor structures, one finds nine independent quark distribution 
functions associating with each tensor structure: 
dX 



— e lXx (PS\^(0)[0, Xw]^(Xw)\PS) 
2n 



q(x,n 2 )p^ + / 4 (x,/i 2 )MV 



(2-18) 



dX 



iXx 



2tt 



(PS\^(0)[0,Xw]r^(Xw)\PS) = 2 Aq(x,fi 2 )p^S-w)+g T (x,fj, 2 )S 



g 3 (x,fi 2 )M 2 w^S-w] 



(2-19) 



dX 
— ( 

2vr 



KX (PS\^(0)[0, Xwjo-^ijs-ipiXw^PS) = 2 [5q(x, n 2 )(Slp" - S^)/M 
+h L {x, fi 2 )M{pV - P V)(5 • w) + h 3 (x, fi 2 )M{S^w u - S u ± w^) 



—e iXx (PS\i}(0){0, \w]il>(\w)\PS) 
2ix 



2 Me(x,fi 2 



(2-20) 
(2-21) 



where we have written 5 M = Sn + S± with Sjt — (S ■ w) p^ + [S ■ p) wA The spin vector 
can only occur linearly in the decomposition, because the dependence on is determined 
by the density matrix (1 + , y 5 $/M)/2 (recall that we normalize so that S 2 = —M 2 ). 

We suppressed the flavor structure in the above definitions. When the quark fields have 
a definite flavor as ip ~ > ipf (f — u,d,s), the distribution functions on the r.h.s. of (2-18) 
- (2-21) should be understood as those for the corresponding flavor : qf(x,fi 2 ), Aqf (x, fi 2 ), 



5qf(x, fi 2 



gl(x,fi 2 ) 



L (x, ), etc. Alternatively, we may insert the flavor matrices A* or 



the unit matrix l in between the quark fields; then the distribution functions on the r.h.s. 
will be understood as those for the non-singlet or singlet part. 
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These quark distributions are dimensionless functions of the Bjorken variable x; they 
depend on the renormalization scale \x as well because the nonlocal light-cone operators on 

the l.h.s. are renormalized at \i. Their spin, 



Table I. Spin, twist and chiral classification of 
the the quark distribution functions. 



Twist 


2 

0(1) 


3 

0(1/Q) 


4 

0(l/Q 2 ) 


spin ave. 


q(x) 


e(x)* 


h(x) 




Aq(x) 


h L (xf 


9s(x) 


s± 


5q(x)* 


9t(x) 





twist and chiral classifications are listed in 
Table I. The distributions in the first row 
are spin-independent, while those in the sec- 
ond and third rows correspond to the lon- 
gitudinally (S||) and transversely (Si) po- 
larized nucleons. Each column refers to the 
twist. The distributions marked with "*" 
are referred to as chiral-odd, because they 
correspond to chirality-violating Dirac ma- 
trix structures r = {0^75, 1}. The other distributions are chiral-even, because of the 
chirality-conserving structures r = {7^, 7^75}. 

In the massless quark limit, chirality is conserved through the propagation of a quark. 
This means that in the DIS one can measure only the chiral-even distributions up to tiny 
quark mass corrections, because the perturbative quark-gluon and quark-photon couplings 
conserve chirality. On the other hand, in the DY and certain other processes, both chiral- 
odd and chiral-even distribution functions can be measured, 0) because the chiralities of the 
quark lines originating in a single nucleon are uncorrelated (see Fig. 2). 

One convenient way to understand the twist classification of distribution functions di- 
rectly from their definitions is to go over to the infinite momentum frame P + ~ Q — > 00, 
so that S ■ w ~ 1 and S± ~ M. This determines the power counting in Q of all terms on 
the r.h.s. of (2-18) - (2-21). The first, second, and third terms in (2-19) and (2-20) behave 
as 0(Q), 0(1) and 0(1/Q), respectively, and the terms in (2-18) and (2-21) are of 0(Q), 
0(1/Q) and of 0(1). 

A conceptually different approach to a similar twist counting is based on the light-cone 
quantization formalism. 0)© In this approach quark fields are decomposed into "good" and 
"bad" components, so that ip = ip+ + ip- with ip + = |7~7 +, an d ip- = where 
7 ± = (7° ±7 3 )/v^2- As discussed in Ref. |l^), a "bad" component ip- introduces one unit of 
twist. Therefore, a quark bilinear i/jTi/j contains twist-2 (if) + rip + ), twist-3 (ip + rtfj_+tfj_rifj + ) 
and twist-4 (ip-Tip^) contributions. This also explains why a unique distribution function 
is defined for each twist and polarization (see Table I): in each of these quark bilinears 
corresponding to twist-2, -3 and -4, the total spin of quark pairs can be or 1. By taking 



Equation (2-17) vanishes for r — 175 due to the time-reversal invariance. 
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the nucleon matrix element, the former describes a spin-averaged nucleon state, while the 
latter describes a longitudinally and transversely polarized state depending on the total 
helicity. The physical meaning of this classification is that a "good" component represents 
an independent degree of freedom. On the other hand, the "bad" components are not 
dynamically independent and can be reexpressed by a coherent quark-gluon pair. 0) Only 
the twist-2 distribution functions are literally the "distributions", which simply count the 
number of independent degrees of freedom having a definite quantum number (flavor, helicity, 
etc.), and correspond to the parton model. In particular, by comparing both sides of (2-18) 
- (2-20), we obtain 

g(x,/i 2 ) = J —e lXx (PS\miOAw}ibij(\w)\PS), (2-22) 
Aq(x,fi 2 ) = j — e iX *(PS\\\m[^MM{M\PS\\)i (2-23) 
5q(x,V 2 ) = ^J ^e i ^(PS ± \m[0,XwW±^w)\PS ± ), (2-24) 

where \PS\\) (\PS±)) denotes the nucleon's spin being in the helicity (S\\ -w = 1) state (in the 
transverse direction). Because the bilinear operator tjjrfjip = (1/ ' P + )i[)^i + ip = (1/ P + )i/)\_i/) + 
corresponds to the number density operator in the light-cone formalism, q(x) gives the 
total number density of a quark with the momentum fraction x of the parent nucleon. 
Due to the additional 75 between the quark fields, Aq(x) is the helicity distribution as 
Aq(x) = q<[(x) — q\{x), where 5 , t(l)( x ) * s the number density of a quark with helicity parallel 
(antiparallel) to the nucleon's spin. If one recalls that S(S) = (l+ r y^/M)/2 is the projection 
operator for a Dirac particle with the polarization S^, one recognizes 5q(x) = q+-{x) —q^(x), 
where q^^(x) is defined similarly to q^n\(x), except that its polarization is defined in 
reference to the transverse direction. We call 5q(x) the "transversity distribution" , following 
Ref. [TE|). On the other hand, the higher twist distributions are the multiparton (quark-gluon) 
correlation functions which contain information beyond the parton model. ©'Ill) This point 
is discussed in detail in §§|5] and []. 

Here it is worth noting some basic properties of the distributions functions: (i) By in- 
serting a complete set of the states J2 n \ n )( n \ between the quark fields of (2-18)-(2-21), and 
using the positivity = (P° + P^)/y2 > for intermediate state |n), we obtain for all 
nine distributions (ft = {q, Aq, 6q, e, hi, gr, /k, #3, ^3}: 

0(x) = , |x| > 1. (2-25) 

(ii) The antiquark distributions q^(x), Aq?(x), h{{x), g f T {x), etc., are defined by (2-18)-(2-21) 
for a flavor / with the nonlocal light-cone operators transformed by the charge conjugation. 
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We find, for <p f = {Aq f , g f T , g£, e f } and x f = W, fl, Sq f , h{, h(}, 



-X) 



-x) 



-X f {x). 



(iii) From the positivity on ql Ax), q[_ __+(x), and the formula q^(x) = q((x)+q[(x) 
qL>(x), we obtain the trivial inequalities < q*(x), \6q*(x)\ < qf(x). 



(2-26) 



2.4. Classification of the gluon distributions 

Our analysis can be extended to the gluon distribution functions. The relevant cut 

diagram is given by Figs. 1 and 2, with the gluon lines connecting the short and long 

distance parts. The gauge-invariant definition of the gluon distribution functions is provided 
by §),§),§) 

2 rd\ 



x 



2tt 



(PS | tr w a G a » (0) [0, Xw] w p G Pv ( A w) \ PS) 

= -lg{x, fi 2 )gT - \aQ{x, fi 2 )ie^ Pa wp(S ■ w) 
- Qst(x, fi 2 )ie^S ±a w p + g A (x, fi 2 )M 2 wV, 



(2-27) 



color matrices, and g± = g 



p^w u 



Table II. Spin and twist classification of the 
the gluon distribution functions. 



where G' MU = d^A" — d u A' M — ig[A^, A v ) is the gluon field strength tensor, "tr" refers to the 

p u w^ is the projector onto the transverse direction. 
The form of l.h.s. is motivated by the fact 
that, in the light-cone quantization with the 
w-A = gauge, w a G afl = {l/P + )d + A». The 
r.h.s. defines four gluon distributions corre- 
sponding to independent tensor structures. 

The spin and twist classifications can be 
inferred directly from (2-27) as in §2.3, and 
are shown in Table II. The distributions 
Q(x), AQ(x),Q 3 t(x) and G±(x) mix through 
renormalization with the flavor singlet parts 
of q(x), Aq(x), gx(x) and f±(x) of Table I, respectively, because of the same twist and spin- 
dependence (see §§3 and ||). On the other hand, there exists no gluon distributions that mix 
with the chiral-odd quark distributions. 

The light-cone quantization formalism again gives some insight into the twist classifica- 
tion in Table II. The gluon fields are decomposed into "good" and "bad" components as 
A^ = Aj_ + A%. The "good" components Aj_ represent the independent degrees of freedom, 
which possess helicity ±1. The "bad" components A 1 ^ can be re-expressed by a coherent 
gluon pair which is coupled to the total spin zero state. This explains the particular pattern 



Twist 


2 

0(1) 


3 

0(1/Q) 


4 

0(1/Q 2 ) 


spin ave. 


g(x) 




Gi{x) 


S\\ 


AQ{x) 






S± 




G3T(x) 
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of the spin dependence for each twist appearing in Table II, and it also implies that the 
twist-2 distributions Q and AQ, which involve only the "good" components, are literally the 
"distributions". From (2-27), we obtain 

2 r dX 

Q[x^ 2 ) = — / — e tXx (PS\trw a G a ^O)[0,Xw}w p Gp fl (Xw)\PS), (2-28) 
x J Ztt 

Ag(x,fi 2 ) = - / ^e iXx (PS\ l \txw a G^(0)[0,Xw)w^Gp^Xw)\PS ll ), (2-29) 



x J 2tt 

where Gp^ = \ep llvp G vp . It is easy to see that Q[x) and AQ(x) represent the total and 
helicity distribution of the gluon, respectively. We introduce the right- and left-handed 
circular polarization vectors = (0, — 1, — i, 0)/v2 and e p L = (0, 1, — i, 0)/\/2, and find 
-G + ^G+ = (G +R YG +R + (G +L ) ] G +L and iG +p G+ = (G +R YG +R - (G +L ) ] G +L . Thus 
G{x) = G](x) + Gi(x) and AQ{x) = Q^(x) — Qi(x), where Gi(i)(x) is the number density 
of the gluon with momentum fraction x and helicity parallel (antiparallel) to the nucleon's 
spin. 

Similarly to the quark distributions, we also obtain: (i) For all four distributions, <fi = 



{G, AQ,QsT,Gi}, the support property ( 2-25 ) holds, (ii) The charge-conjugation transfor- 



mation gives, for even = {AG, Gzt] and <f) odd = {G, £ 4 }, 

0even(-a;) = 0even(^) , 4>odd{-x) = ~(f) odd (x). (2-30) 

(iii) From the positivity of Gi,i(x), we obtain G{x) > \AQ(x)\. 

2.5. The three-particle correlation functions of twist- 3 

Coherent many-particle contents of the nucleon are described by multiparton distribution 
functions. In this paper, we will explicitly deal with the twist-3 quark-gluon correlation 
functions, which are defined as 

^^ e ^^-«)(P5|^(0)y*[o,C«;]^(Cw)[Cw,Ati;]^(A«;)|P^ 



p a e^ Pli S ±p V(x,x',n 2 ) + ---, (2-31) 
Mx+iax'-x) (P5|^( ) 7 « 75 [o, (w]gG^((w) [Cw, Xw]^{Xw)\PS) 
-tp a [p u Sl - p"Sl\§{x, x', fi 2 ) + ■ ■ ■ , (2-32) 



dX d( 



e 



i\x+i((x' —x) 



(PS\m° a %(MgG^((w)l(w,Xwty(Xw)\PS) 



2tt 2tt 

* e rcir :jj_ r r ;;_!_,. j 

+ \v a v v 9^ - p^p'gT - p a p*g± + ^gTW^{x, x',^ 2 )^ 



ie aKp p^gt p - p v gU{S ■ w)M$(x, x', /i 2 



(2-33) 
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where the dots stand for Lorentz structures of the twist higher than 3 (the Dirac matrices 
1 and 275 do not produce any new twist-3 distributions). Similarly to the quark distri- 
bution functions of §2.3, the correlation functions can be defined for each quark flavor /, 
giving (x, x'), \P* (x, x'), etc., or for singlet and non-singlet parts by inserting appropriate 

flavor matrices between the guark fields. 



Table III. Spin and chiral classification of the 
quark-gluon correlation functions at twist- 
3 (0(1/Q)). 



spin ave. 


$(x, 




S \\ 


<P(x. 


x'f 


s± 


W(x,x') 


W(x,x') 



These correlation functions describe interfer- 
ences between the scattering from a coherent 
quark-gluon pair and from a single quark. 

The spin and chiral classifications of the 
twist-3 quark-gluon correlation functions are 
shown in Table III: those marked with "*" 



are chiral-odd, while the others are chiral- 
even. By inserting a complete set of intermediate states between the quark and gluon 
fields, one can obtain the support properties of $ and and their interpretation in the 
partonic language. The variables x, x' and x'—x have the physical meaning of the momentum 
fraction carried by the quark, antiquark and gluon respectively; positive values correspond to 
emission of a parton from the nucleon, and negative ones for absorption of the corresponding 
antiparton (see Ref. |T9|)). The functions and $ vanish unless |x| < 1, \x'\ < 1 and 

\x — x'\ < 1. They also have definite symmetry properties under the interchange x <-> x', 
that follow from the parity transformation combined with the time-reversal transformation. 
For S s = and E A = {^,$}, 

S s (x,x') = S s (x',x) , S A (x,x') = -S A (x',x). (2-34) 

Although these correlation functions assume the rather complicated expressions (2-31)- 
(2-33), they in fact describe the amplitude for finding a triplet of particles (q,q,G) with 
simple physical configuration in the nucleon: the total spin of qqG can be or 1 so that it 
gives a nonzero nucleon matrix element. The spin-1 state of qqG with helicity ±1 is possible 
for the transversely polarized nucleon, and is described by the chiral-even functions \P and 
. Because of the chiral-even structure, q and q have opposite helicities, so that their total 
helicity is with total spin and 1 for & and respectively. Therefore, the total spin 1 of 
qqG is exactly that of the transverse gluon in while both the gg-pair and gluon have spin 
1 in \P . This discussion can be made explicit by considering the combination (see (2-31) and 
(2-32)) 

l - (&(x,x')±§{x,x')) Si = 
_. r d\cK eiXx+iC{x ,_ x) / ps 

J 27T27T \ 



^(0)mCw}gG^(Cw)w u [Cw,Xw}^^^Xw) 



PS 
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(2-35) 



where L denotes the "left-handed component" introduced in §2.4. A similar expression with 
L — > R and 75 — > —75 can also be obtained. On the other hand, it is possible to find two 
other relevant configurations of qqG, the spin-1 state with helicity and the spin-0 state, 
for the longitudinally polarized and unpolarized nucleons. The corresponding amplitudes 
are given by the chiral-odd functions $ and $ , respectively. Due to the chiral-odd structure, 
q and q have the same helicities, so that the gg-pair has the spin 1. In (p, the spin of the 
gg-pair and that of the transverse gluon are antiparallel, and therefore qqG has total spin 
0. In <£>, the spin of the gluon is combined with that of the gg-pair to give total spin 1 and 
helicity 0. These considerations also reveal on a physical basis that the four functions in 
Table III constitute a complete set of the twist-3 quark-gluon correlation functions. 

We note that the quark-gluon correlations <P(x,x'),^>(x,x') and &(x, x'), &(x, x') have 
exactly the same spin, twist, and chiral classification with e(x), Hl{x) and gr{x), respectively 
(compare Tables I and III). As discussed in §§|| and || below, one can actually express 
the twist-3 quark distributions e(x), and gr(x) in terms of a certain integral of the 

corresponding twist-3 quark-gluon correlation functions; only such an integral (over the gluon 
momentum) is potentially measurable in inclusive reactions like the DIS and DY. However, 
we will show that the description of the QCD evolution of e(x), h^x) and gr(x) requires full 
knowledge of the quark-gluon correlations <P(x, x'), <P(x, x') and ^(x, x')^{ functions 
of the two variables x and x'. 

The treatment in this section can be extended to the case of the twist-3 three-gluon 
correlation functions, which are relevant to the (singlet) quark distribution gx{x) and the 
gluon distribution Qst(x). Although the extension is straightforward, the actual analysis 
is complicated due to participation of three identical particles, which have color degrees of 
freedom and obey the Bose statistics. We do not go into the details of the treatment of the 
three-gluon nonlocal light-cone operators,©'© but we discuss equivalent results based on 
the local three-gluon operators in §|5|. 

§3. Relation to the local operator approach and renormalization 

In the previous section, we have defined and classified various parton distribution func- 
tions in terms of the nonlocal light-cone operators. For the structure functions to which 
the OPE can be applied, the above definition, of course, leads to the same results as those 
obtained using the OPE. We briefly discuss this connection employing an example in the 
DIS. Based on this connection, we introduce the renormalization group equations and their 
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solutions, which describe the Q 2 -evolution of the parton distribution functions. We also 
remind the reader of an important issue regarding the factorization scheme dependence of 
the parton distribution functions. 

3.1. Equivalence to OPE and renormalization group equation 

Let us first consider the familiar example of the twist-2 term of the unpolarized structure 
function Fi(x, Q 2 )- To reveal the connection between the nonlocal operator approach and 
the traditional OPE, it is convenient to work in the moment space defined by 



M n (Q 2 ) = f 1 dxx n - 2 F 2 ( Xj Q 2 ) . (3-1) 
Jo 



The factorization formula ( |2-4| ) with i = (f, f, G) and A = iV(nucleon) gives us [ substituting 
ffMZ, V 2 ) = Q f {£, V 2 ), ff/N^, V 2 ) = /" 2 ) and f G/N {t, V 2 ) = G(£, /i 2 ) introduced in §§2.3 
and 2.4 (see also (gT§)] 

M n (Q 2 ) = £ C[ (Q 2 /v 2 , a s ^ 2 )) vltf) + C G n (Q 2 /^ 2 , a s ^ 2 )) v^ 2 ) 

f 

2 C s n (Q 2 //, 2 , a s ^ 2 )) vl^ 2 ) + C G n (Q 2 /^ 2 , a s ^ 2 )) ^(/i 2 ) 



9 

+ 6 



\C« S (QVf 2 , *.(!> 2 )) (»< 3 V) + i"« 8, (^)) , (3-2) 



where we have defined 

10 



vl(fi 2 ) = I dxx n 1 \q f (x,ji 2 ) + g /_ (x,/i 2 )| , f^(/i 2 ) = / dxx n 1 Q(x,ii 2 ) 



o 



(3-3) 

C^ G (Q 2 //i 2 , a s (fi 2 )) = jf 1 dxx n - 2 H LG (x, Q 2 //i 2 , a s (/i 2 )) . (3-4) 

Here f = u,d,s for three flavors, and we have used Hf = Hf, because the DIS is a charge- 
conjugation even process. Note that Hf (£, Q 2 /fi 2 , a s (/i 2 )) = (Q^ 2 )// — 1) + 0(a s )] 
and H G {£,Q 2 /fj, 2 ,a a {n 2 )) = 0(a s ), so that C£ = (Q (e ° 2 )// (1 + 0(a s )), and C G = 0(a s ). 
The second expression of (3-2) corresponds to the singlet (t>^(/z 2 ), v G (/x 2 )) and the nonsinglet 
(v^ 3 ^(/x 2 ), (// 2 )) decomposition (see the last paragraph of §2.2): v\ = v% + v jjj + i>*, t^ 3 ) = 
^n - v ni V n^ = v n JrV n~^ v n^ an d an d C^ s are defined by factoring out the flavor structure 
so that C% NS = l + 0{a s ). 

In the traditional OPE approach, (3-2) is known as the moment sum rule. & The C 3 n 
are called Wison's coefficient functions (short distance parts), and the are given by the 



*' In general, C% ^ C^ S . The coefficient functions corresponding to and Vn^ coincide up to the 
flavor structure because of the flavor independence of the perturbative quark-gluon coupling. 
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nucleoli matrix elements of the local composite operators. In the present case, the relevant 
twist-2 local operators are 



= i n - 1 ^ 7 {w£) w • • • D^ip - traces , (3-5) 
Q ^ n = .„-2 G . {w ^ 2c2 . . . D ^ Cn G^} a _ traces , (3-6) 

where _D M = <9 M — igA^ is the covariant derivative, and { } denotes the symmetrization over 
all Lorentz indices. Here "—traces" represents for the subtraction of the trace terms to make 
the operators traceless, which will be suppressed in the following. of (3-6) is in the adjoint 
representation, and the color indices Cj are explicitly shown. In the OPE, the v n are defined 
by 

(PS\0^\PS) = 2v n P^ ■ ■ -P^ or v n = • ■■w, n (PS\0^\PS) , (3-7) 

where we have used P-w = 1. For 0(^" M ™, corresponding to the flavor structure, becomes 
v^,v n and v^'( 8 \ respectively. Similarly, is given by 

(PSlO^c^lPS) = 2 v^P" 1 ■ ■ ■ P^ . (3-8) 

Now let us calculate of (3-3) using the quark distribution function (2-22) based on the 
nonlocal operator approach. Using the basic property of the distribution functions ( |2-25|) 
and (|2-26|) , we obtain (suppressing the flavor index), 

r+oo 



v n = I dxx 11 l q(x) 

oo 




n-1 



2tt5{\)J (PS\i}(0)yb[0,\w}i;(\w)\PS) 

= i (PS\${0W (iw ■ D) n ~ l m\PS) , (3-9) 

where n = even in the DIS. Comparing this result with (3-5) and (|3-7Q , one realizes that 
our definition of the quark distribution function, (2-22), gives results equivalent to the OPE. 
It is easy to see that of (3-3), combined with (2-28), also gives the matrix element of the 
local composite operator (3-6), and is equivalent to the OPE results. For other structure 
functions, one can show that the situation is the same. In general, the moments of the parton 
distribution functions are given by the matrix elements of the corresponding gauge-invariant, 
local composite operators (see §§5 and 6).S 



*) By extending vl of (3-3) as — J a dxx n ~ 1 {qf (x, /j?) + (— l) n qf(x, /J 2 )], the result (3-9) holds for odd 
n as well. 

**' An exception is the first moment of AQ{x). This will be discussed in §5. 
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The bare nonlocal light-cone operators contain ultraviolet divergences due to loop cor- 
rections in QCD perturbation theory. The renormalization of the ultraviolet divergences 
induces the dependence of the parton distribution functions on the renormalization scale /x, 
which is governed by the renormalization group (RG) equation for the corresponding non- 
local operators. In partonic language, this is the Dokshitser-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equation,© and reads 



[i—q J {x,{i , 

Gt/i 



7T 



y 



d 



a s (n 2 



71 



l dy 

x y 



P Gf> 



X 



) q f ~'(y,v 2 )+P, 



X 



«s0 2 ) G{y^ 2 



X 



>" s (/i 2 ) ) q f '(y^ 2 ) + Pgg ( -,a s (/i 2 ) ) G(y,n 2 ) 



(3-10) 



where the indices / and /' run over quarks and antiquarks of all flavors. The Pij with 
i, j = G,u,u,d,d, ■ ■ ■ are the DGLAP-kernels (splitting functions) given by perturbative 
series in a s (/i 2 ). In the dimensional regularization in 4 — 2e dimensions, these kernels show 
up as the pole residue in e for the loop corrections to the corresponding parton distribution. 
Physically, Py gives the probability density of parton decay process j —*■ i + (anything). 
Thus the equations in (3-10) have the simple physical interpretation as the master equations 
of the probability density for partons. Although they can be obtained straightforwardly by 
renormalizing the corresponding nonlocal light-cone operators,©^ the actual procedure is 
rather complicated. For our purposes, it is more convenient to discuss the equivalent results 
based on the local operator language, utilizing the one-to-one correspondence between the 
nonlocal and local operator approaches established above. 

Taking the moment of (3-10), we go over to the RG equations for the corresponding 
local composite operators. It is convenient to write down these equations for singlet and 
nonsinglet combinations. For the singlet channel, we obtain (a = q, G) 

d 



. 



(3-11) 



b=q,G 



Here, 7^ is given by (an appropriate linear combination of) the n-th moment of the DGLAP- 
kernel Jq 1 dxx n ~ l (a s / tx) P ab (x , a s ), e.g., 



7nKO )) 



GG 



dxx r 



-l« s (/^ 2 



7T 



GG 



(x,a s (/x 2 )) 



(3-12) 
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In the local operator language, 7^ is the anomalous dimension (matrix), which describes the 
mixing between the singlet part of the quark operator (3-5) and the gluon operator (3-6) 
under renormalization. Equation ( |3-11|) can be derived directly from the definition of the 
renormalized composite operators: 



u v,s 

U V,G 




u v,s 

U V,G 



where 0^"^" is the singlet part of the quark operator (3-5), the suffix R (0) denotes renor- 
malized (bare) quantities, and (Z n ) a i, represents the corresponding renormalization constants. 
From the fact that the bare operators do not depend on /1, we immediately obtain (|3-11|) 
with 



d 



ab 



E (z v 

c=q,G 

Similarly, we obtain the RG equation for local nonsinglet operators as 



ac M j ( Z n ) cb 



(3-13) 



where k = 3,8. The anomalous dimension 7^ for the nonsinglet part is also related to the 
n-th moment of the DGLAP-kernels, and it is defined through the corresponding renormal- 
ization constant as 

/i d 



-Z 

Z n d[i 



V vv 



NS 



nW'"' 1 " 
V.NS 



where Oy 1 ^" is the nonsinglet part of (3-5). Note that the RG equations ( |3-11| ) and ( |3-18| ) 
decouple for each n because the Lorentz spin n is a trivial invariant under the evolution. 

The coefficient functions in (3-2), C^(Q 2 //i 2 ), obey the RG equations "conjugate" to 
(|3 • 1 1| ) and ( |3-13|) , respectively, so that M n (Q 2 ) is independent of (arbitrary) scale /i. For 
phenomenological applications, it is convenient to take fi 2 = Q 2 in (3-2). Then the Q 2 - 
dependence of M n (Q 2 ) is dominated by the moment of the parton distribution functions 
(nucleon matrix element of the local operators) renormalized at Q 2 ,© which is given by 
solving ( [3 -111) and (|3-13| ) as 

rQ dK 



<(Q 2 ) 
vl k) (Q 2 ) 



E 

b=q,G 

exp I 



T(Q) exp 
Q dK 



Qo K 



K 



%Mn 2 )) 



7„> s (k 2 )) K fc) (^ 



<(Q 2 o) 



ab 



(3-14) 
(3-15) 



where T(q) denotes the "time ordering" with respect to Q, and Q is the starting scale of the 
perturbative evolution. On the other hand, the coefficient functions C J n (1, a s {Q 2 )) depend 
weakly on Q 2 only through a s {Q 2 ). 
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The coefficient functions C ] n (1, ot s {Q 2 )) of (3-2) as well as the anomalous dimensions 
7^(a< j (/i 2 )) and 7^ \a s {ii 2 )) of (3-14) and (3-15) are calculable using perturbation theory. 
The results can be expressed as power series in a s . We have 



C n(0) + JZ C n(i) + 



where C, 



S,NS 
n(0) 



1' ^n(O) 



/n(0) ^ 47T W?l(1) 

(see the discussion below (3-4)) and 

.V 



(3-16) 



a 



In K) = 7Z 7n(0) 



47T 
47T 



7n {OL s ) = — 7n(0 ) + ( — J 7n(l) + • • • > 

\ 2 

J 7„(i) + ••• • 



The running of a s (/i 2 ) is driven by the /3 function 



djj, 



-fa 



16tt 2 



+ 



(3-17) 



(3-18) 



(16tt 2 ) 2 

where fa = 11 - 2JV//3 and A = 102 - 38.W//3, with N f the number of flavors. When 
we substitute only the leading term of fl3-16|) -( [3-18D into (3-2), (3-14) and (3-15), we obtain 
the leading order (LO) prediction for the (Revolution of M n (Q 2 ) in the RG improved per- 
turbation theory. By including also the next-to-leading order terms (one-loop term of the 
coefficient functions ( |3-16|) , and the two-loop term of the anomalous dimensions (3-17) and 
the (3 function ( p-18|) ), we obtain the next-to-leading order (NLO) prediction for M n {Q 2 ). 
In the NLO approximation, Q3-18| ) is solved to give the running coupling constant, 



a,(Q 2 



g\Q 2 



Air 



4vr A,m(Q 2 /4 CD ) 



1 - 



Aln(ln(Q 2 /4 CD )) 
/3 2 ln(Q 2 /^ CD ) 



(3-19) 



where /Iqcd is the QCD scale parameter at two loops and in the corresponding results of 
(3-14) and (3-15), the leading ([a s ln(Q 2 /Q 2 )] n ) and the next-to-leading (a s [a s ln(Q 2 /Q 2 )] n ) 
logarithms are correctly summed to all orders. For the nonsinglet part (3-15), the integration 
over k is straightforward at the NLO level and gives 



,(*) 



(Q 2 ) = L 7 «(°) 



/2/3o 



a s (Q 2 ) - a s (Q 2 ) fa 



1 + 



7 n (0) 



47T 



,(*) 



(Ql) 



(3-20) 



fa V 2#l 2 fa 

where L = a s (Q 2 ) / a s {Ql) . If we set fa — > anu. j n m 
LO results. The NLO result for the singlet part (3-14) is somewhat complicated due to the 
mixing, but it is not difficult to show that, at LO, (3-14) gives 



and 7% -> 0, (&T% and ( |3~20D reduce to the 



v a n (Q 2 



b=q,G 



L^(o)/ 2/3 ° 



ab 



V b n (Ql) 



(3.2i; 



^Iqcd is renormalization scheme dependent (see §3.2 below). 
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We note that, for M n (Q 2 ) of ( |3-1| ), all quantities necessary for the prediction of the NLO 
Q 2 -evolution are known. 

The correspondence between the nonlocal light-cone and local operators like (3-9) holds 
for other twist-2 distribution functions, the helicity distributions, Aq(x, fi 2 ) and AQ(x,/i 2 ), 
and the transversity distribution, 5q(x,fi 2 ), with appropriate substitutions. The scale- 
dependence of these distributions is also governed by an evolution equation similar to (3-14), 
(3-15), ( |3-2(J| ) and ( |3-21p with the anomalous dimensions for the corresponding local operators 



(see §§5 and 6). 

A similar correspondence between the nonlocal light-cone operators and the local compos- 
ite operators is taken over by the three-parton correlation functions. Manipulations similar 
to (3-9) give, for V(x,x') of (2-31), 

Si ^ dxdx'x^x'^&ix,^) 

= (PS x \$(0)fi(iw ■ D) l - l g& u (<d)w v (iw ■ D) fe ~V(0)|P5 ± ), (3-22) 



with k, I = 1,2,.... Similar results can be obtained for \P(x, x'), x') and <P(x, x') of 
(2-32) and (2-33). Therefore, the double moment of the three-parton correlation functions 
correspond to the gauge-invariant local three-particle operators. In particular, (3-22) and the 
corresponding expression for \P(x,x') generate a set of three-particle operators, which coin- 
cide with those obtained using the OPE for the transverse spin structure function g2(x, Q 2 ) 
(see §§. 

The scale-dependence of the three-particle correlation functions is governed by the RG 
equation for the corresponding three-particle light-cone operators. It gives a generalization 
of the DGLAP equation into the three-body case, which is schematically given by 

H^3(x, x', /i 2 ) = f 1 dy f 1 dy'V(x, x'; y, y'; a s (fi 2 ))S(y, y' , ^ 2 ), (3-23) 



X J X 



for the flavor nonsinglet channel, up to terms proportional to the quark mass. Here E(x, x') 
generically denotes (P(x, x'), <P(x,x') and some linear combination of {^{x, x'), ^(x, x')}, 
and V is the corresponding kernel obtained by renormalizing the ultraviolet divergence 
of the relevant quark-gluon nonlocal operator. As in the twist-2 case, it is convenient 
to work in the moment space, based on (3-22). The moments with different sums k + / 
correspond to different spins, and thus they do not mix with each other. Introducing 
f2 nl (fx 2 ) = / dxdx'x n ~ l ~ 2 x' l ~ 1 S(x, x', /i 2 ), we obtain (/, j = 1,2, ... ,n — 2) 



/i|-^(/i 2 ) + £ [r n (a s ( f i 2 ))] i . f2 nj (/i 2 ) = 0. (3-24) 
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Here r n (a s {ii 2 )) is related to the double moment of the kernel V. In the local operator lan- 
guage, it is the anomalous dimension matrix describing the renormalization mixing between 
the relevant quark-gluon local operators 

Q nl ~ (PS\fy)A{iw ■ D) l - l gG vp w p (iw ■ D) n - l ~ 2 ^\PS) , (3-25) 

with A some Dirac matrix structure. The suffix I here labels many independent local opera- 
tors having the same spin n, and the operators with different I for the same n are allowed to 
mix with each other. Note that the number of independent operators increases with spin n. 
We emphasize that this sophisticated mixing originates from the fact that the three-particle 
correlation functions depend on the two variables x and x', and that it is characteristic of 
higher twist distributions. 

The solution of ( 3-24|) is obtained from (3-14) and ( |3-21| ) with the formal replacement 



v® — > Q ri i and 7„ — > r n . The mixing matrix r n can be obtained by renormalizing the 



corresponding local composite operators of ( 3-25 ) following the standard procedure. In the 



course of the renormalization of higher twist operators, however, we encounter some novel 
phenomena which we discuss in §4. 

3.2. Scheme dependence of distribution functions 

Once the anomalous dimensions (3-17) (or the DGLAP kernels of (3-10)) are known, 
one can predict the Q 2 -dependence of the parton distribution functions with (3-14) and 
(3-15). However, the observable cross sections are the convolution of parton distribution 
functions and the hard scattering parts (short- distance parts) (see (3-2), ( j2-4|) and Q2-6|) ). 
Here an arbitrariness comes into play, beyond the LO, in defining (separating) the parton 
distribution functions and the hard parts. This arbitrariness is called the "factorization 
scheme dependence", and is inherent in the factorization procedure of mass singularities 
(collinear divergences) . & In the local operator language, this corresponds to an arbitrariness 
in defining the renormalized local operators. Both the parton distribution functions and 
the hard parts depend on the scheme which one adopts, and only the convolution of them 
that corresponds to an observable quantity is scheme independent. 

The scheme dependence can be stated in the following way in the example discussed in 
§3.1. Let us concentrate on the singlet part given by the second line of (3-2). Taking a 
different scheme is equivalent to the following change (we take /x 2 = Q 2 ): 

<{Q 2 )^< a {Q 2 ) = E K b «Q 2 ))v b n (Q 2 ) , 

b=q,G 



*' The finite terms of the renormalization constants are in general arbitrary, and are fixed by specifying 
a scheme. 
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C:(l,a s (Q 2 ))^C' n a (l,a s (Q 2 ))= £ C c n (l, a s (Q 2 )) (U^(a s (Q 2 )) 

c=q,G 

Here, U is an arbitrary matrix whose components are functions of a s (Q 2 ), and we set 
C% = C% and = C^. Clearly, the moment M n (Q 2 ) does not change as a result of 
this replacement, but the scale dependence of v'^ is now controlled by the RG equation 
( |3 • 1 1| ) with a "new" anomalous dimension, 



{ll)ah -> (7 V/ i 



lab 



ab 



Thus, the matrix elements (parton distribution functions), as well as the hard parts, can 
be redefined simultaneously from one scheme to another. In view of this, it is important to 
use the same scheme for the parton distribution functions as the corresponding hard parts 
when one predicts an observable quantity. From the phenomeno logical viewpoint, we must 
specify the scheme to be adopted when we derive the parton distribution functions from the 
experimental data. 

It is worth noting here that the parameter Aqcd of (|3-19|) , and equivalently the QCD 
coupling constant, is also a "scheme dependent" quantity. This dependence comes from the 
renormalization (not factorization) scheme in which the coupling constant is defined. In 
contrast to the factorization scheme dependence which cancels between the hard parts and 
the distribution functions, the renormalization scheme dependence of the coupling constant 
remains and is unavoidable in the perturbation theory. The perturbative QCD prediction 
with /Iqcd fixed unambiguously (up to neglected higher order corrections) requires at least 
NLO accuracy. 



]4. Renormalization mixing for gauge invariant operators 



Before proceeding to the detailed study of each polarized structure function and its QCD 
evolution, we discuss the renormalization of the gauge invariant local, composite operators 
from a general point of view. We demonstrate that the renormalization mixing among 
composite operators obeys a particular pattern in QCD (and more generally in non-Abelian 
gauge theories). This knowledge will be useful to derive sophisticated mixing among the 
three-particle operators relevant to the Q 2 -evolution of twist-3 structure functions. 

A local composite operator O is renormalized by determining the corresponding coun- 
terterms C such that arbitrary Green's functions with insertion of O + C, 







T 1 [[ ( j )ai (x i )(0(y)+C(y)) 



i=i 







(4-1) 
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become finite. Here 4> ai {xi) generically denote the (renormalized) elementary fields in the 
Lagrangian, (p ai = {■?/>, A^, ■ ■ •}; the subscript a« indicating various fields as well as Lorentz 
and internal symmetry indices is suppressed for simplicity in the following. We employ the 
minimal subtraction (MS) scheme in the dimensional regularization in 4 — 2e dimensions 
for the discussion in this section. Then C is given by the sum of a power series in 1/e 
with appropriate composite operators as coefficients, and these operators mix with O under 
renormalization. 

As is well known, operators that have common canonical dimension and quantum num- 
bers in general mix with each other. & A particular case of this pattern of "operator mixing" 
is represented by those operators which are invariant under some transformation, and this 
case can be easily understood by the Ward-Takahashi identities: Suppose that we have a 
transformation 8 t that leaves the action S invariant, and that the operator O is invariant 
under the same transformation as 8 t O = 0. A Ward-Takahashi identity for ( |4T|) reads 

n I n 

= £ ( T (5 t <f>(x k )) J] <Kx t ){0 (y) + C(y)) 

k=l \ i^k 

When 8t4> is again an elementary field, e.g. 5t<j)(xk) oc <f)(xk), as in ordinary global transfor- 
mations, the first term of ( |4-2| ) is finite by our construction of the counterterms. Therefore, 
the second term is also finite; but, it is finite only if S t C(y) = 0, since C is given by power 
series in 1/e. This means that only those operators, which are invariant under S t participate 
in C and thus mix with O under renormalization. 

To explore the operator mixing in the renormalization of gauge invariant operators, which 
we denote as Oi (i — 1, 2, • • •), however, the above simple argument should be modified in two 
points (we employ a covariant gauge). First, due to the gauge fixing, the action for quantum 
theory is invariant under the BRST transformation, not under the gauge transformation. 
Therefore, 5 t should represent the BRST transformation #brst i n the present context. Note 
that the gauge invariant operators are BRST invariant, but the converse is not true. The 
BRST invariant operators involving the Faddeev-Popov ghosts are not gauge invariant, and 
such operators are now allowed to mix with gauge invariant operators under renormalization. 
These gauge noninvariant operators include those which vanish by the equations of motion 
(see ( |4-3| ) below). We call the rest of gauge noninvariant operators, which do not vanish 
by the equations of motion, the "alien operators" denoted by Bj (j = 1,2, • • •). From the 
nilpotency of the BRST transformation, <5| RST = 0, the basis of the alien operators can be 
chosen so that they are all BRST exact; & i.e., they can be written as £>j = ^brst^ with Bi 
being some local operator. 

The second modification comes from the fact that the BRST transformation of an ele- 



T[] </>(xi)5 t C(y) 



i=l 



(4-2) 
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mentary field £brst0 is a composite operator. Therefore, it is possible that the first term of 
( |4-2j ) ceases to be finite. If so, we must have #brstC 7^ in order to satisfy (|4-2|) . Fortunately, 
these BRST variant counterterms turn out to be particular operators, as can be understood 
from the following observation: Because the BRST transformed composite operators <5brst0 
are themselves finite operators, the possible divergence in the first term of ( (4-2p occurs only 
when Xk = y, and this divergence should be canceled by the second term. Now we argue 
that such cancellation can be realized when C(y) is given by the composite operators £i{y) 
(i = 1, 2, • ■ ■), which are proportional to the equations of motion: 

SS 



£i(y) = FMy)) 



my) 



(4-3) 



Here Ti is some function of <ft, and from this point we call {£i} the equations-of-motion 
(EOM) operators. It is easy to see that a Green's function containing the operators {Si} 
reads 







i=i 



0\= iJ2S^(x-x k )(0 



i+k 







(4-4) 



by integrating by parts in the functional integral over 0. This indicates that <5brstC(?/) 
can generate the divergent terms for x k = y in the second term of ( |4-2|) when C contains the 
Sj with <5brst^' 7^ 0. Therefore, we conclude that the BRST variant EOM operators, as well 
as the BRST invariant EOM operators, now mix with {Oi}. 

Summarizing, under the renormalization of gauge invariant operators {Oj}, the BRST 
invariant alien operators {Bi} and the EOM operators {Si} mix. tlJ Here {Si} involves both 
BRST invariant and variant EOM operators. The renormalization matrix describing the 
entire mixing among these three kinds of operators obeys a characteristic triangular pattern, 
and it is expressed schematically as 



(4-5) 
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ZqB 


Zoe 
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Zbb 
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The reason for this particular form is simple: Both the bare and the renormalized EOM 
operators should vanish when we use the equations of motion, and thus Zso = Zsb — 0- On 
the other hand, for the alien operators, a Ward-Takahashi identity tells us 







i=i 











T (<5 B rst4 (?/)) 



i=l 







*^ The contact term due to 5Fi(4>(y)) / 6<f>(y) oc S^(0) vanishes in 4 — 2e dimensions. 
**) No other types of operators mix. For a complete proof, see, e.g., Refs. ||) and 
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£ o 

k=l \ 



i^k 







Here, suppose that we determine the counterterms C(Bj) for Bj so that the second line 
becomes finite at y ^ Xi (i — 1, ■ ■ ■ , n). Then, from the first line, we see that ^brst^' + 
C{Bj)) = Bj + <5brstC (Bj) is a finite operator, and thus SbrstC (Bj) gives the counterterms 
for Bj, except for the points y = Xi where {£j} participate. Therefore, Z^o = in the third 
row of the above matrix. 

One important point for actual applications is that the physical (on mass-shell) matrix 
elements of both and {Si} vanish, and these operators do not contribute to the final 
results. The proof is straightforward. Using the LSZ reduction formula, 0* a matrix 
element of E j between "physical" states can be written as 



(P\£i(y)\P') 



lim 

on— shell 



d 4 zd 4 z'e ip - z e- iP '- z ' 



P 2 - M 2 ) (P 12 - M 2 ) (0\TJC\z)£i(y)lC(z')\0) . 



(4-6) 



Here K, denotes an interpolating field to create a hadron state \P) with mass M, or, at the 
level of pure perturbative calculations, JC denotes an elementary field <fi to create a state 
with an on-shell quark or an on-shell (transverse) gluon. With the help of ( [4-4[ ), the Green's 
function on the r.h.s. becomes 

(0|T/Ct(z)^)/C(z')|0> ~ 6^(y - z)r(z, z') + 6^(y - z')f(z, z') , 

where r and f are the two-point functions. Based on their spectral representation,!^ r and 
f cannot have a double pole which cancels the on-shell factors of (4-6). Thus we obtain 



<p|6(y)|p') = o 



(4-7) 



Similarly to the situation with (4-6), an on-shell matrix element of Bi is again proportional to 
the Green's function (OlTSi/Ct/ClO). This can be rewritten using a Ward-Takahashi identity: 



(oim/ct/qo) = (o T (<w&) K)K 











TBiS 



BRST 







(4- 



When /C is a gauge invariant interpolating field, ( |4-8[ ) vanishes. On the other hand, K can 
be an elementary field 0, which gives a BRST transformed composite operator 5brst0 on 
the r.h.s. of (f4-8p. Because #brst0 involves a ghost field, it cannot produce (perturbatively) 
physical particle poles to cancel the on-shell factors of (4-6). Therefore, we obtain generally 



(P|B i (y)|P'> = 0- 



(4-9) 



The disconnected, forward-scattering term (0\£i(y)\0) vanishes in 4 — 2e dimensions. 



28 



How is the above general argument relevant to the structure functions? At the lowest 
twist level, the above complicated operator mixing does not come into play, because there 
exists neither an EOM nor BRST invariant alien operator of twist-2. (The EOM as well as 
the alien operators always have smaller spin by at least one unit than the possible highest 
spin operators of the same dimension (see §§5 and 6).) Therefore, only gauge invariant 
fermionic and gluonic operators mix with each other under renormalization, as in §3.1. 

At the higher twist (> 3) level, both EOM and BRST invariant alien operators participate 
in the renormalization mixing with the gauge invariant operators. This is a characteristic 
feature of the higher twist operators.© Because of ( |4-7| ) and ( |4-9| ), however, only the gauge 
invariant operators contribute to the matrix elements corresponding to the moments of the 
structure functions. Namely, only Zqq in the renormalization matrix ( f4-5p is of physical 
importance. From this one can obtain the anomalous dimensions necessary for the prediction 
of the (^-evolution: 

loo = ^ j- ln (Zoo) ■ 

This would suggest that, by treating only the on-shell matrix elements of composite operators 
to work out renormalization, we would obtain enough information without considering the 
EOM and BRST invariant alien operators. However, the calculation of the on-shell matrix 
elements in terms of purely perturbative Feynman diagrams in massless theory introduces 
another complexity. To these matrix elements, not only the one-particle- irreducible (1PI) 
but also the one-particle-reducible (1PR) diagrams contribute. It often happens, especially 
for the latter class of diagrams, that the infrared singularity coming from the collinear 
configuration cannot be regulated,© so that the perturbative calculation becomes subtle 
and potentially dangerous. © To avoid this risky situation, one of the best ways is to evaluate 
the (off-shell) Green's functions like fl4-l|) . In this case, renormalization can be performed 
by evaluating the 1PI diagrams only, and off-shell external momenta provide an infrared 
cutoff. The calculation is quite straightforward with the manifest Lorentz covariance being 
maintained. One price is that we must take into account mixing of the EOM and BRST 
invariant alien operators; but, with the help of the general results discussed above, we 
can assess unambiguously the necessary renormalization constants. We will discuss explicit 
examples for the twist-3 operators in §§5 and 6. 

§5. Chiral-even structure functions 

In this section, we discuss the chiral-even polarized structure functions of twist-2 and -3, 
which can be measured in the DIS. The QCD evolution is considered with an emphasis on 
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the Bjorken, Ellis- Jaffe sum rules and the twist-3 contribution to g2{x,Q 2 ). 



5.1. Twist expansion 

The chiral-even polarized quark and gluon distribution functions are given by (2-19) 
and (2-27), up to twist-4. These distribution functions have been related to the structure 
functions gi(x, Q 2 ) and g2(x, Q 2 ) by Q2-16Q in the free field theory (for a single flavor case); the 
results can actually be made correct at the LO level by incorporating the LO Q 2 -dependence 
into the distribution functions on the r.h.s. As mentioned in §2.2, the definition of the twist 
based on power counting in 1/Q is not exactly the same as that based on the representation 
of the Lorentz group ("twist = dimension — spin" of the relevant operators). As a result, g% 
(gx) receives contributions from both the twist-3 and twist-2 operators. 

Let us explain how to identify the twist-2 and twist-3 operators© for #2 (fi'r)- (A similar 
argument is possible for the twist-3 gluon distribution Q^t of (2-27).) Taking the moment of 
(2-19) with respect to x and applying manipulations similar to (3-9), 

\w„ 1 ---w» n _ 1 {PS\R a ^--^\PS)= P a {S-w) J^dxx^Aqix^ 2 ) 

+ SI f^dxx n - l g T (x^ 2 ) , (5-1) 
where we have kept only terms through twist-3, we get the local composite operator 

#t{mi~a»»-i} = i'^^Y^D^ 1 ■ ■ ■ D^-^tp . (5-2) 

Here the subtraction of the trace terms g^ il _ lj R a ^ 1 '" fln ~ 1 ^ and ga^R 1 ' >Mn ^ , which generate 
twist > 4 terms, should be understood, and the Lorentz indices /ij are symmetrized (compare 
with (3-5)). However, the index a on the 7 matrix is not symmetrized nor antisymmetrized. 



In order for the operator ( |5-2| ) to have a definite twist (spin), it must be decomposed into 



operators which each has a definite symmetry with respect to the Lorentz indices: 



^>{CT/il---£t„_l} _|_ _ 



n-1 

(-1 \ T-x-T-T/d ---/;„ 1 X 
n 

1=1 



(5-3) 

Here pi "^ 1 '"^™- 1 } an( j R{<?{ni]-~nn-i} are twist-2 and twist-3 operators because they have 
spin n and n — 1, respectively. The nucleon matrix element of these operators can be 
parameterized as, according to their definite symmetry, 

(PS\R {a(11 ^ n - l} \PS) = 2a n S {a P^ ■ • -P^-i> , 

(PS\R lcT{ ^ ] -^- l} \PS) = 2 1 ^—^d n (S ff P^ - S^P a ) P^ ■■■ P^- 1 . (5-4) 

n 
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Substituting these expressions into (5-1) and combining the result with the relation Q2-16 ), 
we obtain (n = odd in the DIS) 



rl If 1 1 

/ dxx n ~ l gi(x) = - dxx n ~ 1 Aq(x) = -a n , (5-5) 
Jo 2 J-i 2 

lr 1 1/1 n — 1 

dxx n ~ l (gi(x) + g 2 (x)) = - dxx n ~ g T {x) = - -a n H d n 

2 J-i 2 \n n 



i . _ n 1 r 1 . _ 1 , , 1/1 n - 1 

in) ■ 

(5-6) 

From the above expression, one can write as (x > 0) 

rl 

92{x) = g 2 w {x) + g 2 (x) = -gi(x) + / — g x {y) + g 2 {x) . 
It is ^2(3^) which receives a contribution from genuine twist-3 operators, 

rl fi 1 

/ dxx^Ux) = — — d n , (5-7) 
Jo 2n 

while the contribution from the twist-2 operators, g^ (x), is called the Wandzura-Wilczek 
part.0 

5.2. The structure function gi(x, Q 2 ) 

We are now in a position to discuss the spin structure function g\. Let us recall the 
moment sum rules for g± based on the local composite operators. The following twist-2 

operators, corresponding to the moments of the nonlocal operators (2-23) and (2-29), are 
relevant here (see also ( |5-2| ) and (5-3)): 

i^-^-i = j»-^ 7 Ky 5 £>w . . • D^- l} ?P , (5-8) 

In these expressions, we employ notation similar to that in (3-5) and (3-6), but the color 
indices are suppressed in (5-9). We define the matrix element of these composite operators 
renormalized at fi 2 by 

(PSlR^f^-'i^lPS) = 2a n (fi 2 )S {a P^ ■ -.P^ , (5-10) 
{PSlRXc'^-'i^lPS} = 2a G ^ 2 )S {a P^ ■ -.P^ , (5-11) 

where a n can be defined for a definite flavor structure (see the discussion below ( |3-7| )). The 
moment sum rules for g\(x, Q 2 ) read (we set /j, 2 = Q 2 ) 

I n (Q 2 ) = f 1 dxx n - l 9l (x,Q 2 ) 
Jo 



1 
2 



2 -E s n (1, a s (Q 2 )X(Q 2 ) + E°(l, a s (Q 2 M(Q 2 ) 



+ ~E™(l,a s (Q 2 )) (af(Q 2 ) + ^)(g 2 ) 



(5-12) 
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where the first (second) line on the r.h.s. gives the flavor singlet (nonsinglet) part, and we 
denote the relevant hard scattering coefficients as E° n , in analogy to the coefficient functions 
of (3-2). Note that (5-12) reduces to (5-5) in the LO for a single flavor case. 

The calculation of the hard scattering coefficients at NLO was performed many years ago 
for both the non-singlet and singlet© parts. The corresponding one- and two- loop© 
anomalous dimensions (DGLAP kernels) have been also calculated. In the following, we 
restrict our discussion to only the first (n = 1) moment, which is related to particularly 
interesting sum rules: the Bjorken sum rule© (the nonsinglet part) and the Ellis- Jaffe sum 
rule© (both the singlet and nonsinglet parts). Equation (5-12) reads for n — 1 (p: proton, 
n: neutron) 



NS, 



l,« s (Q 2 )) 



1 



(8), 



+ ±E?(l,a s (Q 2 ))AZ(Q 2 ) 



with 



(3) 



,(8) 



"1 > 



+ a? - 2a[ , AE = al + <4 + a\ 



a? 



Note that the gauge invariant gluon operator (5-9) (the corresponding coefficient function) 
does not exist for n = 1 in a gauge invariant renormalization scheme, which we assume here 
and in the following. Using Ef s (l,a s ) = Ef(l,a s ) = 1 — a s /n to the one-loop accuracy, 
the NLO expression for If ' n becomes 



12 



(Q 2 



a s {Q* 



71 



(5-13) 



We stress that all systematic QCD analyses of experimental data have been performed at 
this NLO accuracy. 

The fermion bilinear operators which contribute to ( |5T3| ) are the axial vector currents 
ipj^j^ij;. The nonsinglet axial vector current is conserved in (massless) QCD, so that a^ 
and of ( |5-13| ) are independent of the scale Q 2 (the corresponding anomalous dimension 
vanishes to all orders in perturbation theory). Therefore, these values can be fixed using 
information obtained at low energy. Since the Bjorken sum rule is for the difference If — If 
and thus receives a contribution from only the nonsinglet channel, we are led to a definite 
QCD prediction for this sum rule in the sense that the relevant (nonperturbative) matrix 

(3) (3) C 

element a\ is known. We have a\ = = 1.26 from the experimental value of the neutron 
/3-decay. At the NLO of the QCD correction, the Bjorken sum rule reads 



if{Q 2 )-lf{Q 2 



IGa 
QG V 



1 - 



a s (Q 2 



71 



For some higher order calculations, see Ref. |37|) 
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This sum rule has been well established by detailed studies EJ ) '!)'I1 ) 'B ) of the experimental 
data. On the other hand, the situation for the flavor singlet channel is quite different than 
that for the nonsinglet case. At the LO approximation, the singlet axial vector current 

is conserved, so that AS of ( |5-13j ) is also independent of the scale Q 2 . However, there is 

(s) 

no available information at low energy to fix this. If we use the information a\ = 0.58, 
which is obtained from the experimental value of the hyperon /3-decay with the flavor SU(3) 
symmetry relation, and "assume" that we can neglect the strange quark contribution a* = 0, 
we are led to the Ellis- Jaffe sum rule, 



12 1 36 1 



The experimental value for this sum rule is much smaller than the above prediction. This 
suggests the importance of the contributions from the strange quark and/or the higher 
order QCD corrections. When one goes beyond LO, the singlet axial vector current is no 
longer conserved, due to the Adler-Bell-Jackiw anomaly. @^ The corresponding anomalous 
dimension 7J 4 starts at the two-loop order© (compare with (3-17)): 



qq = 12C 2 (R)N f ( 



a, N - 



+ ••• . 

Here C^-R) = 4/3 for N c = 3, and Nf(= 3) is the number of flavors. This makes AS(Q 2 ) 
of ( |5-13| ) scale dependent.©'^ Because there is no operator mixing in this moment, the 
corresponding NLO Q 2 -evolution is obtained by p-20] ) with ^(o) ~^ 0, Tn(i) — > 12C2(R)Nf. 

The scale dependence is weak because it is purely a two-loop effect. 

Concerning the disagreement between the experimental data and the Ellis- Jaffe sum 
rule, and the so-called "spin crisis" problem, namely AS ~ 0.2 from the analysis of the 
DIS data, which corresponds to the total spin carried by the quarks in the naive parton 
model, the role of the polarized gluon has received much attention. As noted above, the 
gauge invariant local operator which would contribute to ( 5T3j) does not exist. On the other 
hand, the polarized gluon distribution AQ(x) has been defined in (2-29) based on a gauge 
invariant nonlocal operator, and their higher moments certainly give the local operator (5-9) 
as / dxx n ~ 1 AQ(x) = for n > 3 (see (5-11)). Where does the first moment of AQ(x) go? 
To resolve this puzzle, let us consider the first moment of AQ(x) in detail. S In this case, the 
situation is somehow subtle, due to the presence of x in the denominator of (2-29). As usual, 
this type of singularity in the distribution functions should be understood as a principal value 



33 



distribution (in the mathematical sense). From the charge conjugation property ( 2-30Q , we 
obtain 

1 

dxAQ(x) 



o 



= I °°dxP(-) [ dXe iXx (PS ll \tiw a G a ^(0)[0,Xw}w p Gf 3lM (Xw)\PS ll } , 

2ll J — oo \X J J 

where P(l/x) denotes the Cauchy principal part. Since J dxP(l/x)e lXx = ine(X), we obtain 
for the r.h.s., 

^ J dXe(X) (P5 ! ||tr WQ G'^(0)[0,AH^ /3 G^(A^)|P 1 S||) . 

This result shows that, in general, the first moment of AQ(x) cannot be expressed by a local 
operator. However, if one chooses the light-cone gauge, w ■ A = 0, we obtain [0, Xw] — > 1 
and w v G v>t (Xw) = (w ■ d)A^(Xw) = (d/dX)A»(Xw), so that, 

£dxAQ(x) = ^JdXe(X)(PS ]l \ti^A^Xw)w /3 G f5tl (0)\PS ll ) 

= - (PS|||tr^(0)^G^(0)|PS||) , (5-14) 

provided that the matrix element (PSyltr A^(Xw)w l3 Gi3 fl (0)\PS\i) vanishes when A — > ±oo. 
Here, tr ^(O)^G^(O) = w^K p /2 with 

Kfs = 2tr^(0) (G^{0) - ^ Act A a (0)A ct (0)) 

being the Chern-Simons current. This argument indicates that the first moment of AQ(x) 
can be expressed by a "local operator" in a specific (light-cone) gauge. Once the first 
moment of AQ(x) is introduced as a finite quantity, it is in principle possible to define a 
schemed where AS of (|5-13|) is replaced by some mixture of a^ + af + al and (5-14), that is 
a{ 7^ AE = a^ + af + af. Among possible schemes of this type, the so-called "AB scheme" 
is often used in the literature. However, it should be stressed that any scheme is acceptable 
on general grounds of QCD, as explained in §4. 

5.3. The structure function g2{x, Q 2 ) 

There have been many works on the transverse spin structure function g2{x, Q 2 ). ©S'&ll!''© 
The main subject in this section regards the derivation of the Q 2 -evolution of its twist-3 part. 
In contrast to the traditional covariant approach based on the local composite operators, 
many authors have employed other approaches to investigate it.B)'©'© These works are 

*) For the phenomenological studies on the g\ structure function and the spin crisis problem, see, e.g., 
Ref. §. 
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based on the renormalization of the nonlocal operators in the light-like axial gauge, or in the 
background field gauge. We develop a framework based on the local composite operators in 
a covariant gauge. 0>'0>il)>li3) According to the discussion in §4, this approach is convenient 
and straightforward. 

To identify the relevant twist-3 operators and a relation among them, one convenient 
method is to utilize the exact operator identities among the gauge invariant nonlocal oper- 
ators. We first recall that we have already obtained a twist-3 local operator R^ivil-p-n-i} Q f 
(5-3), which we denote as R^p Atn_1 : 



•n-1 



■ttn 



n,F 



n 



(n - 1)^7 <T 75-D {m ■ • • D^- l} 4> 



n-1 



1=1 



(5-15) 



As mentioned in §2, there also exist other gauge invariant operators, the number of which 
increases with spin (moment of the structure function), and they mix with each other. 
To infer these additional operators, we note the following operator identity, which can be 
obtained by explicit differentiation (for the derivation, see Refs. [14]) and |51])): 







dz. 



^(O) 7 CT 7 5 [O,.#02) 



^(0)7"7 5 [0,;#(z) 



dt$(0)[0,tz] /|z 75 (t - - ) gG°r(tz)z p - -gG°<\tz)z p \ [tz, z]^{z) 
+ 2m^(0)7 5( 7 CT %[0,z]^) 



+ 



V>(0) 75( x CT %[0, z\{%jp- m q )iP(z) - $(0)(i# + m q )^a ap z p [0, z]^(z) 



(5-16) 



This holds for each quark flavor. Here m q represents the quark mass generically. This 
identity is exact through twist-3. We have neglected the twist-4 contributions and the total 
derivatives which are irrelevant for the forward matrix elements. 

Taylor expanding (5-16) around z p = 0, it is easy to see that the l.h.s. generates operators 
(5-15) with n = 2,3, ■ • ■. Similarly, we find the following twist-3 operators from the r.h.s. 
(flavor matrices and subtraction of trace terms suppressed): 

1 



K n,l 



L n,m 



2-,(TUl---^l n _l 

n n,E 



, W-K-i-z + ty + t^-i-O (Z = l,---,n-2) 
-i n Sm^^D^ ■ . . D ftn -*'f n - 1 ')J) , 



2n 

+ ^7 CT 7s^ Atl 



■4>{i IP - m (? )7 <T 7 5J D Ml • • • £>m»-2 7 m«-i<0 



(5-17) 
(5-18) 

(5-19) 
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where 



y x = -i n g Si)D^ ■ ■ ■ G aH ■ ■ ■ D^Y"- 1 ^ , 
Ui = -i n - l gS4>D^ ■ ■ ■ G af1 ' ■ ■ ■ D^- 2 Y n '^ ■ 

Here, S denotes the symmetrization over «j. The R n ^ correspond to the first term of (5-16); 
it can be shown that this term is related to the integral of a quark-gluon correlation function 
E(x, x') = &(x, x') + W(x', x) + &(x, x') - W(x', x) (see (2-31) and (2-32)), and thus the R n j 
are given by the double moment of S(x, x') in analogy to (3-22). The R n . : E are the gauge 
invariant EOM operators, and the R n , m represent the quark mass effects, which cause the 
mixing of the chiral-odd operators. In the flavor nonsinglet case, neither the pure gluonic 
operators nor the BRST invariant alien operators is involved, and (5-15) and (5T7)-(5-19) 
form a (over) complete set of the twist-3 operators (except for the BRST noninvariant EOM 
operators, as discussed below). 

From the above Taylor expansion of (5-16), we recognize that the local operators (5-15) 
and (5-17)-(5-19) satisfy the relation Eil 

71-1 ™" 2 

KT^ 1 = — JCT""- 1 + E(» - 1 - + KT^ 1 ■ ( 5 ' 2 °) 

n i=i 

This equation implies that these operators are not all independent. Therefore we have 
freedom in the choice of the independent operator basis to work out renormalization. We 
stress here that we can take any basis, and the final results do not depend on the choice. 
A convenient choice of the independent operators is (5-17), (5-18) and (5-19). For the n-th 
moment, n gauge- invariant operators participate in the renormalization for the nonsinglet 
channel. 

The RG equation for the matrix elements of the local composite operators read (compare 
with ( gg§ ) 



where 



fi-ffnk(^) + Y,[ r nM^))] k . iW(/i 2 ) = , (k,j = !,•••, n-2, m) (5-21) 



j 



(PSlRlf'^i^lPS) = 2f nl (fx 2 ) (S^P" 1 - S^P a )P^ ■ ■ -P^- 1 , 
(PS\RZ^" tln - 1 (^)\PS) = 2/„ m (u 2 ) (S a P^ - S^P a )P^ ■ ■ -P^- 1 . 



It is possible to derive a relation like (6-17) and ( |6-18| ) below among Qt^x), 5(x,x / ), m q Sq(x) and 
Aq(x) from (5-16). 

lis operator identity can also be obtained by applying the relation [D^, D v ] = —iqG^ v to 

M 



*) 

(5-15) 
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Equation (|5-21 ) at the LO is solved to give, similarly to ( |3-21 



fnk{Q 2 ) =J2[ L 



(5-22) 



where we defined the mixing matrix X n as F A (a s ) = (a s /2ii)X n . 

The anomalous dimension (matrix) i - ^ (ct s (/i 2 )) can be obtained by calculating the renor- 
malization constants for the n operators ((5T7)-(5-19)) according to the general argument 
in §4. The calculation proceeds in the standard way. We multiply the operators by the 
light-like vector to symmetrize the Lorentz indices and eliminate the trace terms. We 
calculate the quark-quark-gluon off-shell three-point Green's functions with the insertion of 
the relevant operators. As discussed in §4, the gauge (BRST) noninvariant EOM operators 
are also required to complete the renormalization. These operators are obtained 0-*'© by re- 
placing some of the covariant derivatives _D Mi by the ordinary derivatives <9 Mi or in (5-19). 
This causes the actual calculation to become very complicated because there are (2 n ~ 2 — 1) 
operators of this kind. This technical problem can be overcome by introducing the vector 
n^, which satisfies w^Q^ = 0,© and by contracting the Lorentz index of the external gluon 
with this vector. This has two merits: First, the tree vertices of the gauge invariant and 
noninvariant EOM operators coincide. Thus essentially only one EOM operator is involved 
in the mixing. Second, the structure of the vertices for the twist-3 operators are extremely 
simplified, and the computation becomes more tractable. An explicit one loop calculation 
appears in Refs. [I7|) and The results are consistent with the general argument in §4, 
namely the renormalization constants take the following form: 
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Rn,l 


\ 




/ 




Rn,m 
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Rn,E 


J 





V 












{Z n )lE 





J 



( 



R 



'I -J 



Rn,m 
Rn,E 



nllm 



The necessary anomalous dimensions are unambiguously determined from (Z n )ij , {2,. 
and (Z n ) mrn . For details, we refer the reader to Refs. (47p and f|8]). 

Now let us consider the moment sum rule for g2(x,Q 2 ). We define the matrix elements 
of composite operator (5-15) by (5-4) with rWM-i^-i} _> i?^"'* 1 "" 1 (fj, 2 ) and d n -> d nF (fx 2 ). 
The moment sum rule for the genuine twist-3 part of g2(x,Q 2 ) reads (we consider a single 
flavor case for simplicity) 

n — 1 



dxx n g 2 (x,Q 



In 



-d nF (Q 2 )Ei(l,a s (Q 2 )) 



*> For the description of other processes like DY, r£ for even n, as well as for odd n, are required (see 
(6-2)), and these anomalous dimensions have recently been worked out in connection with renormalization 
of light-cone wave functions. l3 
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1 

+ 2 



/.m(Q 2 W(l,a,(Q 2 )) + £ UQ 2 )E' n (l,a,(Q 2 )) 



n-2 

(5-23) 

1=1 



where E*(l, a s {Q 2 )) is the hard scattering part corresponding to each operator. Note that, 
at LO, (5-23) is consistent with ( |5-7| ) (see below). From the constraint Q5-20| ), the matrix 
elements satisfy 



n-1 . .... n-1 



-d nF (Q 2 ) = f nm (Q 2 ) + J2(n-l- l)f nl (Q 2 



n n l=1 



It is to be noted that the explicit form of the Q 2 -evolution of each term in (5-23) de- 
pends on the choice of operator basis. The anomalous dimensions as well as the hard parts 
E*(l, oi s {Q 2 )) take different forms, depending on the basis. Let us demonstrate this point at 
LO: in the basis of independent operators, which includes R n<F , it turns out that the hard 
parts become 

Ef(l,0) = l, E™(1,0) = ^(1,0) = 0, 
and the result coincides with (|5-7|) . On the other hand, if we eliminate R n ,F, we have 

— Ti — 1 — 

^(1,0) = — , £*(l,0) = n-J-l. 

However, the moment (5-23) itself remains the same, as should be the case. 

Now the present discussion reveals that the (^-evolution of the moment (5-23) obeys a 
quite sophisticated renomalization mixing due to ( |5-22j ). Even if we neglect the quark mass 
operator (5-18), (n — 2) independent terms contribute to the evolution of the n-th moment, 
and mix with each other. As a result, the comparison of this (Revolution with experimental 
data is very difficult; the numerical diagonalization of the (n — 2) x (n — 2) mixing matrix, as 
well as estimates of the (n — 2) independent nonperturbative parameters. However, we note 
that in the iV c — > oo limit, this complicated renormalization mixing is drastically simplified: 
Each component of the mixing matrix of (|5-22 ), [X n ]^ (k,l = 1, • • • , n — 2), involves two 



Casimir operators, C 2 (R) = N £ N 1 and C 2 (A) = N c . In the large N c limit, that is, with the 
replacement (^(i?) — > 4r, ^ ne following exact relation holds: H* 

n-2 

£( n _ k - 1) [X n ] kl = (n-l- l)w n , (5-24) 

k=l 

where w n is the smallest eigenvalue of X n in this limit and is given by 

( n ~ 1 1 1 1 \ 
m n = 2N c + — . (5-25) 
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As a result, the moment obeys simple Q 2 -evolution without mixing. 

f 1 dxx n -%{x, Q 2 ) = L mn/I3 ° ^ dxx n - x g 2 {x, Q 
Jo Jo 



(5-26) 



Therefore, in this limit, the twist-3 nonsinglet structure function g 2 (x,Q 2 ) obeys the Q 2 - 
evolution of the twist-2 nonsinglet type, and involves one single nonperturbative parameter. 
We discuss further implications of this N c — > oo simplification in §6.2, in connection with 
similar simplifications of the chiral-odd twist-3 distribution functions. 

In the flavor singlet case, the situation becomes much more complicated. We have opera- 
tors that are made of two or three gluon fields, as well as the BRST-invariant alien operators, 
in addition to the above set of operators bilinear in quark fields (with the flavor matrix being 
removed) : 



1 n,G 

rp<Jil\---il n -\ 

1 n,l 



1 n,B 



1 n,E 



i n_1 5 A S 
i n ~ 2 g S 



G V ^D^ ■■■D^GJ 7 

n — 1 

1 = 2 — — 



(5-27) 
(5-28) 



•n-l 



S{G afl1 D^ 2 ■ ■ ■ £)^"- 2 } a 



1 



a 



•n-l 



S{G U ^D^ 



+ gf ahc {d^t)x c ) 
D^- 2 } a {(D v G v ^- 1 ) a 



+ #i a 7 ^-^ + -d^-HSTAl) - gf^id^-^X }- 
a 



(5-29) 



(5-30) 



Here A antisymmetrizes a and fi±. The color structure of these operators should be under- 
stood similarly to (3-6). \ and x are the ghost fields, t a is the color matrix in the fundamental 
representation, [t a , t b ] = if abc t c , and a is the gauge parameter. T nj ; is trilinear in the gluon 
field strength G^ v and the dual tensor, and thus represents the effect of the three gluon cor- 
relations. It satisfies the symmetry relation j^ 1 '"^™- 1 — T^^j Mn_1 due to the Bose statistics 
of the gluon. T n ^ is the gluon EOM operator which is BRST invariant. T n b is the BRST 
invariant alien operator, which is the BRST variation of the operator, 



•n-l 



S{G a>11 D^ 2 ■ ■ ■ D^ n - 2 ^ a Q^ n - 1 ^ a 
The gluon bilinear operator (5-27) is related to the trilinear ones (5-28) by§)>§ 

\l-2)Cr 2 {n - I - 2)C n n Zt 



I n,G 



n-l 
2 

E 

1=2 



2(n-l) 2(n- 

n{n - 2l)C?~? 
2(n-2)(n-l) 



+ (-1) 



i+i 



1 V rp<Jfli---fl n -i 
~ L ) I n,l 
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i=l 



(5-31) 

where Cy = n!/[r!(n — r)!]. To derive this relation, we have used [D^, D u ] = ~igG^ v and 

a^x vol ~t~ D Q: a ~\~ D a G uu 0, H) a G va ~\~ D U G aa ~\~ D Q G au ^vaapD \G 

where the first identity is the usual Bianchi identity and the second one is a consequence of 
the relation 

£a/3v8 + 9p8 

As a result of (5-31), Ell we can conveniently choose a set of independent operators as (5-28) 
- (5-30). For the n-th moment, these independent operators mix with each other and 
with the n gauge-invariant operators bilinear in the quark fields discussed above under the 



renormalization. till An explicit calculation is given for the n = 3 case in Ref. [19|), and the 
result is again consistent with the general argument in §4. 

Finally, let us comment on the Burkhardt-Cottingham sum ruleS which corresponds to 
the first moment of Q2(x,Q 2 ). Due to the lack of local operators for n — 1 moment, the 
operator analysis suggests the sum rule 



/ dxg 2 (x, Q 2 ) = . 
Jo 



The validity of this sum rule has been a long-standing problem.® However, it was shown 
that this sum rule is satisfied in QCD perturbation theory, lit 



§6. Chiral-odd structure functions 

This section is devoted to discussion of chiral-odd distribution functions of twist-2 and 
-3. As mentioned in §2.3, these distributions are inaccessible to totally inclusive DIS because 
of their chiral property, and thus there exist no experimental data for them. These "new" 
distributions are expected to open a new window to explore the nucleon spin structure, 
and several experiments designed to access them have been proposed. For example, the 
twist-2 distribution 5q(x, /J, 2 ) of (2-20) and (2-24) can be measured through the polarized DY 



When one takes the physical matrix element of (5-31), the EOM and the BRST invariant alien oper- 



ators drop out ((4-7) and (|4-9|)) , and (5-31) becomes the identity suggested in Ref. 



**' Up to the gauge (BRST) noninvariant EOM operators. 
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process with two transversely polarized proton beams, where the LO double transverse-spin 
asymmetry Att reads 



£/(Q (e ° 2 )// 


'5qf(x l ,Q 2 )Sqf(x 2 ,Q 2 ) + (f< 




E/(Q^ 2 )// 


g'(*i,W>2,Q 2 ) + (/^ 


►/)" 





Att — a TT r — > (6"1 



up to 0(M 2 /Q 2 ) corrections. Here / = u,d,s,---, citt is the partonic asymmetry corre- 
sponding to the short- distance part, and we set fi 2 = Q 2 with Q the dilepton mass. The 
variables X\ and x 2 refer to the parton's light-cone momentum fraction in the two nucle- 
ons labeled as "1" and "2", respectively. Other processes to measure 5q(x,Q 2 ) include the 
direct photon production, two-jet production, and heavy-quark production in transversely 
polarized nucleon-nucleon collisions, il)>Il) and also semi-inclusive DIS observing a certain 
hadronic final state with one0> or two pionsi^ or a polarized A® in the current fragmen- 
tation region. 

On the other hand, the twist-3 distribution functions /^(x, /i 2 ) and e(x,fi 2 ) of (2-20) and 
(2-21) usually constitute small corrections to the leading twist-2 term, so it is difficult to 
extract them experimentally. For instance, e(x,fi 2 ) contributes as an 0(M 2 /Q 2 ) correction 
to ( |2-6| ) for the DY process between two unpolarized nucleons.l^ However, 1il(x, /i 2 ) is some- 
what immune to this difficulty, like the chiral-even twist-3 distribution gx{x, /i 2 ) discussed in 
the last section. h L (x, fi 2 ) reveals itself as a leading contribution to the polarized DY process 
between the longitudinally and transversely polarized proton beams. The corresponding LO 
double spin asymmetry Alt is given by© 

A LT = 

E f (Q {el)2 )ff [Agf(x 1 ,Q 2 )x 2 g T (x 2 ,Q 2 )+x 1 h{(x 1 ,Q 2 )Sgf(x 2 ,Q 2 ) + (/ ~ /)" 
E/(Q^ 2 )// W> 2 ,Q 2 ) + (/ <- /)] 

(6-2) 

up to corrections of twist-4 and higher. Here the nucleon "1" ("2") is longitudinally (trans- 
versely) polarized, and a L T is the partonic asymmetry. Actually, o>lt is of 0(M/Q), while 
a TT = O ((M/Q)°). Thus Alt is smaller by a factor of M/Q compared to Att- This is sim- 
ilar to the situation with the chiral-even structure function g 2 (x, Q 2 ) compared to gi(x, Q 2 ). 

These chiral-odd distribution functions will be measured at BNL, DESY, CERN, etc., 
in the near future.© In particular, RHIC at BNL will provide us with the first data of 
Sq(x, Q 2 ) and Hl{x, Q 2 ) through Att and Alt- In view of this, it is essential to work out a 
detailed study of these distribution functions based on QCD, including Q 2 -evolution. © 
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6.1. The transversity distribution 5q(x,fi 2 ) 

The transversity distribution Sq(x,fi 2 ) constitutes a complete set of the nucleon's twist- 
2 quark distribtion functions together with the spin-independent distribution q(x, fi 2 ) and 
the helicity distribution Aq(x, fx 2 ) (see (2-22)-(2-24) and Table I). As demonstrated in §2.3, 
Sq(x, fi 2 ) has a probabilistic interpretation as Sq(x, fi 2 ) = q<-(x, fi 2 )—q^(x, fi 2 ) corresponding 
to a simple parton model. For nonrelativistic quarks, where boosts and spatial rotations 
commute, one would obtain <5g(x,/i 2 ) = Aq(x,fi 2 ) (compare (2-23) and (2-24)). However, it 
is well known that the quarks inside a nucleon cannot be nonrelativistic, so that 8q{x, /i 2 ) 
and Aq(x,fi 2 ) should contain different and complementary information with regard to the 
nucleon spin structure. 

In the local operator formalism, the transversity distribution is described by the ma- 
trix elements of a set of twist-2 composite operators, which are obtained by following the 
procedure outlined in §3.1. Taking the moment of (2-20) and using ( 2-25Q , we obtain 



/l r+oo 
dxx n ~ l 5q(x, fx 2 ) = / dxx n ~ 1 b~q(x,fi 2 ) = t n (fi 2 ) 
-1 J — oo 



(6-3) 



where n = 1, 2, 3, 



and 



(PS\0 



T 



(fi 2 )\PS) = —t n {fi 2 )S {S V P^ - S^P V ) P» 2 --- P^ 



with <y 



(fi 2 ) the twist-2 local operator renormalized at /i 2 : 

i^Sija^i^D^ 2 ■ ■ ■ D^tp 



(6-4) 



(6-5) 



Note that the subtraction of all the trace terms is implied, and t n can be defined for a definite 
flavor structure as t{ (/ = u, d, s) or P n (j = q, (3), (8)) (see the discussion below (|3-7|)) on 
the r.h.s. of ( |6-4| ) and ( |6-5| ). In particular, the first (n = 1) moment t\ is called the "tensor 
charge" . For a quark of flavor / and in the rest frame of the nucleon, (|6-3] )- (p-5| ) give 



dx [5q f (x, /i 2 ) - 5q f {x, fx 2 )) = t{{fi 2 ) , {PS\^ } a^f\PS) = 2t{{fi 2 )S j 



where we have used ( |2-26P and set <jj = crjoi^. It is instructive to compare (|6 
corresponding expression for the axial charge a\ (see (5-5) and (5-10)): 

J dx (Aq f (x,fi 2 ) + Aq f \x,fi 2 )) = a{(fi 2 ) , (PS\^ f a^ f \PS) = 2a{{fi 2 )S j 



(6-6) 
with the 

(6-7) 



Here, the relevant operator corresponds to the spin operator for a Dirac particle, and differs 
by 7° from that of (|6-6| ). Therefore, t\ is certainly different from a\ for relativistic quarks. 



Compared with the notation of (3-5) and (5-8), we suppress the subscript "ip" in (6-5), because the 



gluonic operator O t G 



does not exist in the present case (see the discussion below (|6-8 



)• 
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Comparison between ( |6-6| ) and ( |6-7| ) also reveals that t\ probes the contribution from the 
valence quarks only, while a\ includes the helicity of the sea quarks. This reflects the 
difference between the charge-conjugation properties, expressed by ( |2-26| ), for 5q(x,fi 2 ) and 
Aq(x,fi 2 ). This discussion suggests that measurements of t 1 , combined with those of a 1; 
will provide us with deeper knowledge about the effects of relativistic and sea-quarks on the 
nucleon spin structure. 

The /independence of 5q(x, n 2 ) is governed by the RG equation for the moments (|6-3|) 
as (compare with Q3-13 )) 



^n(^) + 7nM^))tn(V 2 ) = , (6-8) 

where 7^ (a s (/i 2 )) is the anomalous dimension of the operator O^' 1 '"''" of ( |6-5|) , and is related 
to the corresponding DGLAP-kernel. As discussed in §2.4, there is no gluon analogue for 
Sq(x,fi 2 ), due to the chiral-odd nature (compare Tables I and II). Therefore, 5q(x,fi 2 ) does 
not mix with gluons under renormalization, in contrast to the chiral-even distributions dis- 
cussed in §5; the evolution of each moment is governed homogeneously by a single anomalous 
dimension j^(a s (fi 2 )), as ( |6-8|) , both for the nonsinglet (t^ 3 \t^) and singlet (t g ) channels. 
We expand 7^(a s ) similarly to (3-17) with one- and two-loop coefficients 7^(0) anc ^ 7n(iv The 
quantity 7^ was calculated H)'© to be 




7^0) = 2C 2 (R) l + 4£- , (6-9) 



and the LO Q 2 -evolution of t n {Q 2 ) is given by ( |3-2(J| ) with 7^ — > 7^ ), ln(i) ~^ 0- ^ n or der 
to demonstrate unusual evolution properties due to (|6-9| ), we compare it with the well known 
result for the one- loop contribution to 7^ of ( |3-11| ), 



v 



7n(0) — 7n(0) 



-y 




2C 2 (R)[1- . ■ rr + 4V^ . (6-10) 



7n(o) wnon y governs the LO evolution of the nonsinglet part of q(x, Q 2 ) as well as of 

Aq(x,Q 2 ); i.e. 7^ (0) = 7^(o)-^ We see that 7«(o) > 7«(o) ^ nolds for a11 n - This id- 
eates that 5q(x, Q 2 ) displays a stronger Q 2 -dependence than the chiral-even distributions 
q(x,Q 2 ) and Aq(x,Q 2 ). In particular, 7^ > while jYto) = 0- Therefore, the tensor 
charge t\ is a scale-dependent quantity, unlike the vector and axial charges; its absolute 
value decreases with increasing energy scale as 

tl (Q 2 ) = L c ^'H x {Ql) . (6-11) 



43 



This reflects the fact that the local operator ipUv^i^^ip relevant to the tensor charge is not a 
conserved current. [ On the other hand, as discussed in §5.2, 7j / (a s ( / u 2 )) = 7^(a s (/i 2 )) = 
to all orders, because of conservation of vector (nonsinglet axial vector) current.] We also 
note that the difference between 7^ and 7^ is larger for smaller n. This suggests that the 
evolution of Sq(x, Q 2 ) has a rather different behavior in the small-x region compared with 
that of q(x, Q 2 ) and Aq(x, Q 2 ), because the small-n region in the moment space corresponds 
to the small-x region in the Bjorken-x space. 

Recently, the two-loop anomalous dimension J^ci) has been calculated©' in the 

MS scheme. It was shown that all the evolution properties of 5q(x, Q 2 ) discussed above 
are preserved and become even more pronounced by including the NLO effects.© It is also 
worth noting the following point to demonstrate the importance of the NLO corrections: the 
rightmost singularity of the anomalous dimension on the real axis in the complex n plane is 
known to determine the small-x behavior of the corresponding parton distribution function 
within the DGLAP evolution (see, e.g., the second reference of Ref. [7])). The rightmost 
singularity of ^(o) * s located at n — — 1, because 2" =1 1/j = ip(n + 1) + j E , where ip(z) 
is the di-gamma function, and 7^ is the Euler constant; this gives 5q(x, Q 2 ) ~ x (ignoring 
logarithms) for x — > 0. On the other hand, the corresponding singularity of Jnfi) is at n = 
Therefore, only after including the NLO effects, the DGLAP asymptotic behavior 
gives 8q(x,fi 2 ) ~ const, and becomes consistent with the Regge asymptotics discussed in 
Ref. [71J. 

Now, the NLO evolution of all the twist-2 distribution functions of the nucleon is avail- 
able both for the nonsinglet and singlet parts. It is possible to perform the QCD analysis for 
physical quantities and to make a systematic comparison between those distribution func- 
tions at the NLO level. As emphasized in §3.2, this NLO evolution must be combined with 
the NLO short- distance parts calculated with the same scheme as the former, to compute 
an observable quantity. For the present chiral-odd case, the corresponding NLO correction 
to ott of ( |6-1|) appears in Refs. [67]) and |7l]), and applied to give a prediction of A^t-© 



Theoretical estimates of the x-dependence of 5q(x,fi 2 ), or equivalently, the matrix ele- 
ments ( |6-4| ), require elaborate calculations based on nonperturbative methods. In this con- 
nection, it is worth noting that 5q(x, /i 2 ) obeys some inequalities: the first one is \5qf(x, /i 2 ) | < 
qf(x,fi 2 ), due to the positivity of the probability density, as discussed below ( |2-26|) in §2.3. 
The second one was proposed by Soffer,^ based on the positivity properties of (diagonal- 
ized) helicity amplitude for the forward quark-nucleon scattering, 

2\8q f (x,v 2 )\ < q f (x,fi 2 ) + A ? / (x,/i 2 ) . (6-12) 

Both inequalities hold for each flavor of quark and antiquark. These inequalities could 
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provide useful constraints on 5q(x,[i 2 ), combined with the known information of q(x, /x 2 ), 
Aq(x } /j, 2 ), and/or with reliable assumptions. We note that there have been several dis- 
cussions in which it has been pointed out that the inequality (|6T2| ) might suffer from higher 
order corrections beyond the LO.B'^'S 

There have been various theoretical estimates of the tensor charge ti based on lattice 
QCD,0 QCD sum rules,© nucleon models,® etc. There also exist calculations of the 
x-dependence of Sq(x,fi 2 ) based on QCD sum rules, and nucleon models. We do 

not go into the details of these results here, but simply note that the results suggest that 
5q(x,fi 2 ) is not small, but is of the same order as Aq(x,fi 2 ). 

6.2. The twist- 3 distributions /i 2 ) and e(x,/i 2 ) 

We now proceed to a systematic study of the chiral-odd twist-3 distribution functions 
h^x,/! 2 ) and e(x, /i 2 ) of (2-20) and (2-21). Our first task is to demonstrate that these distri- 
bution functions contain information regarding the quark-gluon correlations in the nucleon, 
and thus that they are quantities beyond the parton model. This can be visualized by re- 
expressing hi(x, fi 2 ) and e(x, (i 2 ) in terms of the quark-gluon correlation functions introduced 
in §2.5. As in the chiral-even case g<i in §5.3, the existence of such nontrivial relations reflects 
the fact that the basis of the twist-3 distribution functions introduced in §2 is overcomplete: 
on account of the QCD equations of motion, the number of independent dynamical degrees 
of freedom is less than the number of distribution functions corresponding to independent 
Lorentz structures. 

The standard technique to reveal such constraints is to derive relations between local 
operators which arise by taking the moments of the definitions (2-20), (2-21), and (2-33), 
as worked out for g 2 (x,Q 2 ) in §§5.1 and 5.3. For hi{x,fi 2 ), a similar approach requires 
somewhat complicated analysis^ of the trace terms of twist-3, which were subtracted out 
from the local operator fl6-5|) . Here we employ another approach, which is more convenient 
for the present case. S We directly derive the relations among the relevant distribution 
functions in the Bjorken-rc space by using exact operator identities between the nonlocal 
operators. The relevant operator identities are (compare with (5-16)) 

{^(o)z 7 5<v^[M^)} 



Soffer's inequality (6-12) cannot be converted into that among the matrix elements t n ,a n , and v n of 



(p-4[), (5-10) and (3-7) corresponding to the moments of the distributions because of the mismatch in the 



charge-conjugation property ( 2-26| ) among 8q(x, /i 2 ), Aq{x, /i 2 ) and q(x,[i 2 ). In particular, it is useless in 



deriving a model-independent constraint on the tensor charge t\. 

**) A similar technique was applied to the analysis of the twist-3 light-cone wave functions of vector 
mesons. E$ 
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= -i J dv (v - tp(0)i-f 5 cr a(S z (3 [0,vz]gG au (vz)z l/ [vz,z}4j( y z) 
+ im q ip(0) ^7 5 [0, z)tp(z) 

+ % - {^(0)(z - m q )fr B [0, z}^(z) + $(0# 76 [0, z\{iip- m q )^{z)} , 

#)) [0,2^) -#))#)) 

(iv , ?/'(0)(T Q/3 z ( g[0, vz]^G Q!l/ (f z)z u [v z, uz]ip(uz) 
—im q ^(ti) / [0, wz]^(ttz) 

- | {^(0)(z p-m q )jt [0, uz]V>(u*) + ^(0) j, [0, ^> - m q )i/;(uz 



(6-13) 











/o 


2 7o 



(6-14) 



These relations are exact up to operators containing total derivatives which are irrelevant 
for the parton distribution functions. It is straightforward to prove the identities directly 
for the nonlocal operators (see also Ref. |51])), or to demonstrate that their formal 
Taylor expansions at small quark- ant iquark separations generate towers of identities for the 
corresponding local operators, which are derived in Ref. [16]) (see also Refs. |53[) and p6|)). 

In the light-cone limit z 2 —>■ 0, the nucleon matrix elements of the nonlocal operators 
on both sides of (6-13) and (6-14) can be expressed in terms of the distribution functions 
defined in (2-18)-(2-21), and (2-33); note that the matrix elements of the operators involving 
the equations of motion vanish due to (|4-7|) . For the calculation of the l.h.s. of (6-13), we 
need the tensor decomposition of the matrix element of the nonlocal operator, where the 
quark-antiquark separation is not restricted to being light-like: 



(PSmO^a^^z^lPS) 
2 r/„ S-z \ 1 



Mir p-z v 



P M ) (P ■ z) I d.n 



-ixP 



- M- 



, S ■ z 
P-z 



^—P 
P-z \ 



dxe 



-ixP-. 



-1 



5q(x) + 0(z 2 
h L (x) - 5q(x) + 0(z 2 ) 



(6-15) 



where the coefficient of each tensor is determined by matching with the light-cone limit 
(2-20). The matrix element of (6-13) now yields the differential equation 



-x 2 -T- (~h L (x) 
dx \x 



25 g (x)-P I" dx> 1 ( 9 d ^^^ m « d 



7l^-^]^0-^^), (6-16) 



-l x — x' \ dx dx 



where "P" denotes the principal value, and we have used the symmetry relation <P(x, x') = 
—(£>(x',x) of (|2-34j) . The solution of this equation with the boundary condition Q2-25| ) reads 



h L (x,fi 2 ) 
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2x 



£{X) dy Sq ( y,fI ^ 



M 



V 

X 



e(x) 



dy-P 



dy — 
i y 



- 2x 



y 



(6-17) 



where — 1 < x < 1 and e(x) 
calculation, 

xe(x, jj 2 ) = P 



X \x\ 



We also obtain from (6-14), by straightforward 



dx' 

- : a; — x 



1 7 ^(x,x',/i 2 ) + m9( ^ ) ^- 



M 



-q[x,fj, 



(6-18) 



where ^(x, x') = <&{x' ', x) of ( |2-34|) has been used. The results (6-17) and (|6-18|) were obtained 
by a different method employing the light-like axial gauge in Ref. |87|). Here we have derived 
them with manifest gauge and Lorentz covariance being maintained. 

The results (6-17) and (|6-18|) show that hi(x, /i 2 ) and e(x, fi 2 ) can be completely expressed 
in terms of the other distribution functions. According to the various terms on the r.h.s. of 
(6-17) and (|6-18|) , we decompose the solution in an obvious way as 

h L (x, /i 2 ) = hf w (x, /i 2 ) + h2(x, /i 2 ) + h£{x, /i 2 ) , 



e(x,/i 2 ) 



e c (x,/i 2 ) + e r < 



(6-19) 
(6-20) 



Here, h^ w denotes the first term of (6-17), which is the contribution from the twist-2 
distribution Sq(x,fi 2 ) and corresponds to the Wandzura-Wilczek contribution for the case 
of g2(x,Q 2 ) (see §5.1). On the other hand, the other terms stand for the effects due to 
the genuine twist-3 operators: hi and e c represent the "dynamical" twist-3 contributions 
expressed as the particular integral of the quark-gluon correlation functions $(x, x') and 
<P(x,x'), whose explicit forms are given by appropriate projection of (2-33) as 

-1 



$(x, x') 



2M 



d^ d( i\x+i£( x '- x ) 
2tt 2tt 



x (P5||^(0) 75 a a %[0,CH^G'^(C^)^[C^,AH^(Aw)| J P5||) , (6-21) 

and similar expression for @(x, x') with 75 — > 1 and S\\ — > S. Furthermore, the quark mass 
m q generates another type of twist-3 effect, h™ and e m , which are given by the chiral-even 
twist-2 operators multiplied by m q . It is straightforward to show that (6-17) and Q6-18Q give 
the following relations© in the moment space (compare with (5-6), (5-23)): 



M n [h L (fi 2 ] 
M n \e(fi 2 ] 



t n {ix 2 )+M n h%(fi 



n + 1 



+ 



n 



n + 1 M 



l n— 1 



Mr 



e<V) 



+ 



M 



-v n -i(fi 2 ) 



(6-22) 
(6-23) 



*' Actually, e(x, /i 2 ) contains the terms proportional to the delta function S(x). We eliminate these terms 
by multiplying e(x,/i 2 ) by x, because they are irrelevant to the partonic interpretation. 
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where t n , a n _i and v n -i are defined as (|6-4j ) , (5-10) and ( |3-7| ), respectively, and 



M n 






f ^dxx n - x h L 


X, /I ) , 


M n 


h c L (n 2 ) 




^dxx n ' l h c L 


(x,/i 2 ) = 


M n 






l^dxx n - 1 e c ( 


x,/i 2 ) = 



-1 



dxx n e(x, /i ) 



21 



1=2 



1=2 



n + 1, 
„ 1-f-l 



(6-24) 



4 



- & n. K +(/^ 2 ) ■ 



Here 



and 







n — 1* 








.2. 




2 



(6-25) 
(6-26) 



n-l-l/l-2 



x n - 2 $(x,x',ti) 



n-l-l/l-2 



x n 2 ${x,x' ,n 2y 



(6-27) 



Substituting (6-21) and the corresponding expression for ^(x, x') into (6-27), 6 ni ; and b n j are 
given by 



(P5|H^j'" M "- 1 |P5') = 2b n:l MS(S >11 P^ ■ .-P^\ 
(PSIW^'^-'IPS) = 2b n jM(P til P» 2 ■ - -P^- 1 ) 



(6-28) 
(6-29) 



with the twist-3 local quark-gluon operators, 



■n— 1 



n.l 



-Sipa a ^~f 5 D^ ■ ■ ■ gG^ 1 ■ ■ ■ D^- 1 ^ -(l^n-l + 1) 



■n— 1 



-Sipa^D^ ■ ■ ■ gG^ 1 ■ ■ ■ D^ip + (I -+ n - I + 1) . 



(6-30) 
(6-31) 



Pl—Pn-l 



On the r.h.s. of (6-28)-(6-31), subtraction of the trace terms is implied. Note that W£\ 
and Mn_1 are odd and even under the replacement I — > n — I + 1, respectively, cor- 
responding to the symmetry of <P(x,x') and <P(x,x') under the interchange x <-> x' . The 
W'5'" Mn_1 with Z = 2, • • • , «~ (W£t'"' <n ~ 1 with Z = 2, • • • , «+) form an independent basis of 
the twist-3 quark-gluon operators for h L {x,ji 2 ) (e(x, /i 2 )) in the local operator approach. 

Next, we proceed to the (Revolution of h^x.Q 2 ) and e(x, Q 2 ) (see (6-19) and (6-20)). 
The Q 2 -dependence of Sq(x, Q 2 ) has been discussed in §6.1, and this completely determines 
that of hY w (x, Q 2 ). Similarly, h™(x,Q 2 ) and e m (x,Q 2 ) are driven by the (Revolution of 
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Aq(x, Q 2 ) and q(x,Q 2 ) discussed in §§5.2 and 3.1, combined with the running quark mass 
m q (Q 2 ). On the other hand, the Q 2 -evolution of the quark-gluon correlation contributions 



h ( l(x,Q 2 ) and e c (x,Q 2 ) is complicated: $(y,y\jj 2 ) and <P(x,x' ', /i 2 ) of (6-17) and (|6-18|) are 
respectively governed by a RG equation similar to ( |3-23| ) with the corresponding kernel [up 
to the mixing of m q (fi 2 )Aq(x, /i 2 ) and m g (/i 2 )g(x, /x 2 )]. 

Following the discussion in §3.1, let us go over to the moment space to treat the evolution 
of the quark-gluon correlations. Neglecting the quark masses for simplicity, {b n> i{ix 2 )} and 
{b n ,i(^ 2 )} of (6-28) and (6-29) obey RG equations similar to (|3-24|) with the corresponding 
anomalous dimension matrices. In the LO evolution, we denote these anomalous dimension 
matrices as f^Y~ and f^Y^ ', respectively, where the superscripts — and + refer to the 
"parity" of the relevant operators (6-30) and (6-31) under / — > n — 1 + 1. Ell The RG equatk 
are solved to give (compare with ( |5-22j )) 



cions 



KAQ 2 ) =J2[ L 



k=2 



Ik 



bn,k{Qc 



bnAQ 2 



k=2 



Y+/f3 



Ik 



b n ,k(Qo) 



(6-32) 



The mixing matrices have been obtained by computing the one-loop corrections to the 
relevant operators (6-30) and (6-31) in Refs. |53|) and|5B|) (see also Ref. S In accord with 
the general argument in §4, the EOM operators, which are generated by Taylor expanding 
the nonlocal EOM operators appearing in (6-13) and (6-14), play roles in the course of 
renormalization. ciiJ On the other hand, in contrast to the chiral-even twist-3 case of §5.3, 
one does not encounter the mixing of the pure gluonic operators nor the BRST invariant 
alien operators in the singlet channel, due to the chiral-odd property. 



and M. , 



Substituting ( ggg ) into the r.h.s. of (6-24) and (6-25), M n [h c L (Q 2 
exhibit a complicated mixing pattern characteristic of the higher twist operators. They in- 
volve k~ — 1 and k+ — 1 independent terms corresponding to the eigenvalues of the mixing 
matrices Y~ and Y+, respectively. More and more terms with different anomalous dimen- 

for n > 6 
by the 



sions contribute to M n h^Q 2 ) (M n e c (Q 2 ) J for larger n, and M n h1(Q 2 ) 
M n e c (Q 2 ) for n > 5) are not directly related to M n h^(Ql) ^M n e c (Q 



evolution, in contrast to the twist- 2 case of §6.1. It would not be possible to distinguish 
experimentally between terms with different anomalous dimensions in Ai n h^iQ 2 ) and 
M n \e c (Q 2 



*' The RG equations at LO are not affected by the insertion of 75 into the relevant operators (6-30) and 
(6-31). @ Therefore, we label the mixing matrices by the "parity" of the corresponding operators. 
**) In these works, the mixing of m g a„_i (m q v n -i) with b n ^ (6 n ,fc) is also calculated. 
***) xhe BRST noninvariant EOM operators also come into play and can be treated© similarly to 
the case of 55.3. 
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There exist, however, two important limits, N c — > oo and n — > oo, where .M n h1(Q 2 
and .M n g c {Q 2 ) obey simple Q 2 -evolution. The mixing matrices of ( |6-32j ) involve two 
Casimir operators C 2 {R) = and C^A) = Af c . In the large N c limit with the replacement 
C 2 {R) 



the exact relations 



1=2 



E 

1=2 



Ik 



1 - 



21 



n 



y: 



2k 



Ik 



a; 



Y~ 



n + 1 

1 - (-1)" 
4 



(6-33) 
(6-34) 



have been derived, E^'Ea) where A^ are the smallest eigenvalue of the mixing matrices in 
this limit, and are given by 



(6-35) 





As a consequence of these relations, we obtain (compare with (|5-26|) ) 

M n \h c L (Q 2 )} = L x -^M n \h c L (Ql)} , M n \e c (Q 2 )} = L xi ^M n \e c (Q 2 ) 



(6-36) 



which indicate that M. n h^lQ 2 



and Ai r 



obey simple Q 2 -evolution without any 



[e C (Q 2 ] 

complicated operator mixing. They are governed by the anomalous dimensions given in 
analytic form ( |6-35| ). The inverse Mellin transformations of (|6-36|) show that /i£(x, Q 2 ) and 
e c (x, Q 2 ) obey DGLAP-type homogeneous evolution equations, as in the twist-2 nonsinglet 
case. 

The mathematical reason for this simplification is the same as that for the similar simpli- 
fication applied for the chiral-even (nonsinglet) structure function (72(2, Q 2 ) in (|5-24|) -( |5-26| ): 
The coefficients of the relevant quark-gluon contributions b n j_ and 6 n) ; in (6-24) and (6-25) give 
the left eigenvector of the mixing matrices Y~ and Y^ corresponding to the lowest eigenvalue 
A^ and A+ (see (6-33) and (6-34)). These phenomena can be referred to as the decoupling of 
the particular combination of the quark-gluon operators, which give /i£(x, /1 2 ) and e c (x,fi 2 ) 
as (6-24) and (6-25), from the evolution ( |6-32| ). The same decoupling is observed at large 
n for arbitrary values of A^.i^'S In this case, we obtain ( |5-24j ), (6-33) and (6-34) with the 
lowest anomalous dimensions ( |5-25| ) and ( |6-35[ ) shifted according to (g n = zu n , A~, A+) 



Qn ^ Q n + (4C 2 (R)-2N, 




(6-37) 



to O (ln(n)/n) accuracy. With these modifications of the anomalous dimensions, the results 
( |5-26| ) and ( |6-36| ) are valid to O (l/N 2 ■ ln(n)/ra) accuracy. 
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We recognize that the simplifications of the (Revolution in the two limits, N c — > oo and 
n — > oo, are universal phenomena for all twist-3 distribution functions g2(x,Q 2 ), hi(x,Q 2 ) 
and e(x,Q 2 ). These results provide a powerful framework both in confronting experi- 
mental data and for model building, since the results are valid up to corrections ~ V-^c' 
which are numerically small. From a general point of view, they are interesting in provid- 
ing us with an example of an interacting three-particle system in which one can find the 
exact energy of the lowest state.® For phenomenology, the description of each moment of 
the twist-3 distributions now requires one single nonperturbative parameter. This is note- 
worthy because it means that inclusive measurements of twist-3 distributions are complete 
[to O (1/N 2 ■ ln(n)/n) accuracy] in the same sense as those of the twist-2 distributions; i.e., 
knowledge of the distribution at one value of Qq is enough to predict its value at arbitrary Q 2 
in the spirit of DGLAP evolution equation. For model building, the results are used ESS 
to rescale the model calculations at a low scale to the large values of Q 2 which correspond 
to actual experiments. It is also possible to predict i># the behavior of the twist-3 distribu- 
tions in the limits x — > and x — *■ 1: The rightmost singularity of the relevant anomalous 
dimensions ( |5-25 ) and (|6-35|) are located at n = on the real axis in the complex n plane, 



and thus the N c — > oo DGLAP evolution equations predict the following behavior for x — > 
with (a s (g 2 )/vr)ln(l/x) < I: 

g^ s (x,Q 2 ), h c L (x,Q 2 ), e c (x,Q 2 ) - const. , (6-38) 

where g2 S (x,Q 2 ) denotes the nonsinglet part of the genuine twist-3 piece of g2(x,Q 2 ), and 
"const." denotes a Q 2 -dependent normalization factor of each distribution. On the other 
hand, the n — > oo DGLAP evolution equations for these distributions are driven by a common 
anomalous dimension ( p-37 ), from which we obtain for x — > 1 with {a s {Q 2 ) / tt) In (1/(1 — x)) 



< 1 

g? s (x,Q 2 ), h c L (x,Q 2 ), e c (x,Q 2 ) * F{^ 2 ) p ( ^ 2) - - - L ^(l-x)^ K , (6-39) 

where K = — (AC-iiR) / fio) InL, r(x) is the gamma function, and je is the Euler constant. 
F{fi 2 ) and c(fi 2 ) are /^-dependent constants, which are specific to each distribution. The 
results (|6-38| ) and ( |6-39| ) give a guide to the expected small x and large x behavior, which is 



important for experimental extrapolations. 



*) Similar simplifications are observed in twist-3 nonsinglet fragmentation functionsEj) and in twist-3 
meson wave functions e3 (see also Ref. p0|)). 
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7. Summary 



In this article, we have studied and surveyed the polarized structure functions in QCD. 
Our starting point was the facrotization theorems in QCD, which provide strong machinery to 
separate systematically a variety of high-energy cross sections into short- and long-distance 
contributions. To describe a "universal" long-distance part, we have defined the parton 
distribution functions as the nucleon matrix element of gauge invariant nonlocal operators, 
and discussed complete classification of quark and gluon distribution functions, as well as 
those of twist-3 three-particle correlation functions. One of the key points of our approach is 
that we preserved maximal (gauge (BRST) and Lorentz) symmetries of the theory at every 
step of investigation. This allowed us to derive the properties and relations of the distribution 
functions in an economical way. Another point is use of the one-to-one correspondence 
between the nonlocal (coordinate or Bjorken-x space) and the local operator (moment space) 
approaches. This again corresponds to the maximum use of the symmetry to simplify the 
dynamics, because the moment n is invariant under the evolution. In particular, we discussed 
in detail the (Revolution of the twist-3 distribution functions and the corresponding three- 
particle correlation functions in the local operator approach, where all the relevant steps can 
be worked out based on the standard and familiar field theory techniques. A characteristic 
feature of higher twist gauge invariant operators was discussed from a general viewpoint, 
and it was shown how we should deal with the renormalization of the complete set of twist- 
3 distribution functions. We also emphasized the essential role of the QCD equations of 
motion to reveal the interrelation between the different twist-3 distributions and to derive 
their (Revolution. 

For all the twist-2 distribution functions, we now have sufficient theoretical information to 
make a NLO prediction of the QCD evolution. Although the treatment of twist-3 distribution 
functions was very complicated already at LO, the corresponding (Revolution has been 
completely worked out giving the mixing among an independent set of nonperturbative 
matrix elements. These results will provide us with a powerful framework to analyze new 
experimental data with sufficient accuracy. As explained in the text, the nature of the parton 
distribution function depends on the scheme employed. Therefore, we need the short-distance 
parts (hard scattering coefficients) in the same scheme for each independent hard process 
to perform a consistent QCD analysis. We need a wide variety of processes to extract a 
consistent set of information and to get an unified picture of the parton distributions. The 
polarized Drell-Yan and other processes utilizing hadron colliders will give us opportunities 
to fully understand the behavior of quarks and gluons inside the nucleon. 

The transversity distribution is the "new" and final twist-2 quark distribution and is ex- 
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pected to play a distinguished role in spin physics, complementing the information regarding 
the nucleon spin structure from the conventional helicity distribution. The "measurable" 
twist-3 distribution functions are also of special interest, because they contain information 
concerning quark-gluon correlations inside the nucleon and are free from renormalon ambi- 
guities. However, the phenomenology of twist-3 distributions is not as straightforward as 
that of twist-2 distributions due to the fact that the twist-3 distributions are essentially 
three-particle correlation functions: We need to determine many parameters to make precise 
predictions for the observable structure functions, so that it is generally very hard to obtain 
full information of the twist-3 distributions. In this respect, drastic and universal simplifi- 
cation of all twist-3 distributions in the iV c — > oo (or n — > oo) limit is crucial, and solves the 
above problem for all practical purposes. Now the N c — > oo and n — > oo evolution equations 
allow for the development of phenomenology and model building of the twist-3 distributions 
in analogy to the twist-2 case to 0(1/N^) accuracy. 

We mention that there still remain many subtleties and controversial aspects which were 
not covered in this article. Even in the twist-2 case, we are sometimes forced to make 
several assumptions on various unknown corrections in using QCD predictions to describe 
experimental data. Among those problems are: the small x (and large x) behavior of the 
structure functions, 0) theoretical treatments of the "wee" partons, including the role of mas- 
sive quarks, the flavor separation of the structure functions, etc. Also, the power corrections 
in 1/Q 2 could produce rather large effects in realistic experiments, and cause some ambi- 
guities. These power corrections come from the higher twist terms as well as target mass 
effects. The treatment of the latter effects is well established,© but, as for the former, it is 
known that there appear ambiguities in the twist > 4 terms due to the renormalon singular- 
ities. All these problems require careful theoretical investigations. Forthcoming precision 
measurements will also provide us with information and clues to answer these questions. 

Finally, we hope that various kinds of new experiments and more detailed theoretical 
investigations will be able to clarify not only perturbative but also nonperturbative aspects 
of QCD related to hadron spin physics. 
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